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ABSTRACT 


TEMPERATURE  MEASUREMENTS  AND  INTERNAL  WAVES 
IN  SENECA  LAKE,  NEW  YORK 

Myron  Harvey  Fliegel 

Temperatures  were  measured  near  the  center  of  Seneca 
Lake  with  thermistors  and  bathythermographs  during  four  dif¬ 
ferent  years.  Seneca  Lake  is  a  long,  narrow,  deep  lake  in 
the  Finger  Lakes  region  of  western  central  New  York  state, 
with  a  temperature  structure  similar  to  that  of  many  oceans. 
The  bathythermograph  records,  which  were  used  to  reconstruct 
the  annual  temperature  cycle  of  the  lake,  indicate  that  it  is 
freely  circulating  throughout  much  of  the  winter  although  its 
temperature  is  below  4°C.  The  annual  heat  budget  was  esti¬ 
mated  to  range  between  40,000  and  47,000  cal. /cm. ^  for  most 
years . 

Thermistor  chains  with  probes  at  five  depths  revealed 
large  amplitude  internal  waves.  Wave  amplitudes  of  25  feet 
were  commonly  seen  and  50  foot  amplitudes  were  occasionally 
seen.  These  internal  waves  are  characterized  by  a  general 
decrease  in  temperature  at  several  levels  for  tens  of  hours 
followed  by  a  very  rapid  rise  that  restores  temperatures  to 
earlier  values  within  a  few  minutes,  followed  by  oscillations 
of  periods  of  tens  of  minutes  lasting  for  several  hours  there¬ 
after.  The  form  of  these  waves  is  very  similar  to  that  of 
non-linear  surface  undular  surges  investigated  theoretically 
using  the  Korteweg-de  Vries  equation  and  with  tank  models  by 
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various  authors.  The  internal  waves  seen  in  Lake  Seneca  are 
therefore  considered  to  be  non-linear  internal  undular  surges. 

A  direction  of  travel  was  determined  for  12  surges  seen 
in  1971  when  two  thermistor  chains,  190  feet  apart,  were  in 
operation.  In  all  12  cases,  the  disturbance  was  travelling 
north.  Five  of  the  surges  were  seen  at  a  third  chain  2.4  miles 
away  and  enabled  phase  velocities,  which  ranged  from  19  to  50 
cm. /sec.,  to  be  determined.  Wind  speed  and  direction  data  were 
collected  for  78  days  in  1970  and  1971.  No  clear  relationship 
between  wind  patterns  or  atmospheric  pressure  patterns  and 
surges  could  be  determined.  It  is  concluded  that  the  surges 
are  generated  by  a  resonant  interaction  between  the  winds  or 
pressure  systems  and  weak  internal  seiches  already  travelling 
through  the  lake. 
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INTRODUCTION 


Introductory  remarks 

During  the  summer  of  1968,  the  Arctic  group  at  Lamont- 
Doherty  Geological  Observatory  installed  a  string  of  five  ther¬ 
mistors  connected  to  a  chart  recorder  near  the  center  of  Seneca 
Lake.  The  initial  purpose  of  the  installation  was  to  test  ins¬ 
trumentation  to  be  used  in  the  Arctic  Ocean  on  ice  island  T-3. 
Lake  Seneca  is  a  deep  fresh-water  lake  with  strong  temperature 
stratification  during  the  summer  that  models  that  of  many  oceans. 
In  addition,  the  Naval  Underwater  Systems  Center,  Systems  Meas¬ 
urement  Branch  at  Dresden,  New  York,  maintained  a  barge  moored 
near  the  center  of  the  lake  which  provided  an  excellent  platform 
from  which  to  conduct  experiments.  For  these  reasons,  the 
excellent  logistics  provided  by  the  NUSC,  SMB  and  the  short  dis¬ 
tance  from  Lamont  (less  than  five  hours  away  by  automobile)  it 
was  felt  that  instruments  could  be  field-tested  at  this  location 
rather  inexpensively  before  being  sent  to  the  Arctic.  However, 
it  became  apparent  soon  after  the  thermistors  were  installed 
in  the  lake,  that  the  temperature  data  obtained  were  very  unusual 
and  extremely  interesting  in  themselves.  Large  amplitude  tem¬ 
perature  fluctuations,  with  periods  of  tens  of  minutes,  were 
seen  at  various  times  on  the  original  records.  When  the  records 
were  digitized  at  one-hour  intervals  and  then  plotted  it  became 
evident  that  the  lake  was  in  continual  motion. 

Temperature  measurements  were  continued  and  expanded  during 
the  following  three  years.  By  the  fall  of  1971,  17  thermistors 
on  three  chains  at  two  separate  barge  sites  were  in  operation. 
This  report  presents  the  results  obtained  from  the  four  years 
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of  data  collection  on  the  lake. 

Regional  setting  and  water  characteristics 

The  Finger  Lakes  of  western  central  New  York  State  lie  in 
an  area  of  gentle  relief  underlain  by  beds  of  limestone,  shale 
and  sandstone  with  a  total  thickness  of  8000  feet  and  dipping 
gently  to  the  south.  These  rocks  are  Paleozoic,  the  surface 
beds  in  the  region  being  Devonian.  Physiographically ,  Seneca 
Lake  lies  principally  in  the  New  York  Limestone  Belt,  with  only 
the  southern  end  extending  into  the  Northern  Appalachian  Plateau. 

The  deep  and  elongated  basins  of  the  Finger  Lakes  were 
formed  from  previously  existing  river  basins  by  glacial  ero¬ 
sion  during  the  Pleistocene  Epoch.  Von  Engeln  (1961)  separates 
the  region  into  east  and  west  provinces  according  to  the  pre¬ 
glacial  drainage  patterns.  In  the  west  the  rivers  flowed  south 
and  the  glaciers  followed  the  valleys  downs lope  in  a  manner 
analogous  to  that  which  created  the  fiords  of  Norway.  In  the 
eastern  province  of  which  Seneca  is  the  western-most  lake,  the 
preglacial  rivers  flowed  north.  Glaciers,  entering  from  the 
north,  had  to  flow  up  the  river  valleys  which  narrowed  to  the 
south.  This  funnelling  effect  enhanced  the  erosive  powers  of 
glaciers  in  the  eastern  province  with  the  result  that  these 
lakes  are  deeper  than  those  in  the  western  province,  Seneca 
and  Cayuga  having  bottoms  below  sea  level. 

Seneca  Lake,  with  a  surface  area  of  67.7  square  miles  and 
a  volume  of  15  x  109  m3  is  the  largest  of  the  Finger  Lakes.  It 
is  35  miles  long  with  a  maximum  width  of  3.25  miles  and  a 
drainage  area  (including  that  of  Keuka  Lake  which  drains  into 
Seneca  Lake)  of  707  square  miles.  The  maximum  depth  is  618 
feet,  which  is  174  feet  below  sea  level,  and  the  mean  depth  is 
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290  feet.  The  lake  surface  is  444  feet  above  sea  level. 

Drainage  is  north  into  Lake  Ontario  by  a  devious  route  through 
the  Seneca  and  Oswego  Rivers. 

Chemically,  Seneca  Lake  belongs  to  the  Northern  Appala¬ 
chian  Plateau  because  its  major  inlets  flow  primarily  through 
that  region  (Berg,  1963) .  Both  Seneca  and  Cayuga  Lakes  have 
much  higher  concentrations  of  sodium  and  chloride  than  the 
other  Finger  Lakes.  Berg  (1963)  considers  the  possibility 
that  this  results  from  discharges  by  salt  mining  operations  on 
the  shores  of  these  lakes  and  dismisses  this  hypothesis  on  the 
basis  of  the  volume  of  salt  involved.  Salt  strata  underlie 
all  of  the  Finger  Lakes  region  with  a  few  of  the  purest  and 
most  extensive  being  commercially  extracted.  The  Seneca  and 
Cayuga  lake  basins  intersect  beds  lower  than  any  of  the  other 
Finger  Lakes.  Berg  (1963)  contends  that  impure  salt  beds  in 
contact  with  the  lakes  are  responsible  for  the  abnormally  high 
sodium  and  chloride  concentrations. 

A  temperature  profile  taken  with  a  bathythermograph  near 
the  southern  end  of  the  lake  (J.  Gorman,  personal  communication) 
showed  a  temperature  maximum  at  a  depth  of  70  feet  that  was  4°F 
greater  than  the  surface  temperature.  This  implies  a  high  con¬ 
centration  of  salt  at  that  depth.  It  was  not  determined  whether 
this  was  due  to  leaching  of  salt  from  a  bed  or  was  due  to  dis¬ 
charges  from  a  nearby  commercial  salt  plant. 

An  account  of  the  geology  of  the  Finger  Lakes  region  can 
be  found  in  Von  Engeln  (1961)  and  an  account  of  the  general 
limnology  of  the  region  in  Berg  (1963) . 


4 


DATA  COLLECTION  AND  HANDLING 
Temperature  probes 

Temperatures  were  measured  using  Fenwal  glass  probe 
thermistors.  Thermistors  possess  a  negative  temperature 
coefficient  of  resistivity,  i.e.,  increased  temperatures  pro¬ 
duce  decreased  resistance.  The  equation  relating  temperature 
and  resistance  may  be  expressed  in  the  form: 

0/T 

R  =  Ae  (1) 

with  A  and  3  constants  for  a  particular  thermistor,  R  the 
resistance  in  ohms  and  T  the  temperature  in  degrees  absolute. 

The  probe's  leads  were  connected  to  the  male  leads  of  an 
underwater  connector  and  then  sealed  in  either  a  brass  or 
plastic  cylinder,  using  a  potting  compound  in  such  a  way  that 
at  one  end  of  the  cylinder  the  temperature-sensitive  part  of 
the  probe  was  exposed.  A  water-tight  seal  was  made  by  connec¬ 
ting  this  to  the  female  connector  and  then  through  cable  to 
the  recording  instrument.  The  time  constant  of  the  thermistors 
used  was  25  seconds. 

The  thermistors  used  had  a  nominal  resistance  of  2000  ohms 
at  a  temperature  of  25°C  and  5800  ohms  at  0°C.  Each  thermis¬ 
tor  was  calibrated  individually  by  accurately  measuring  the 
resistance  of  the  potted  thermistor  and  its  temperature  while 
immersed  in  an  isothermal  bath.  The  temperature  standard  was 
a  Dymec  quartz  crystal  thermometer  which  had  been  calibrated 
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against  a  triple-point  bath.  Values  were  determined  for  at 
least  12  different  temperatures  for  each  probe.  A  least-square 
method  was  used  to  determine  the  constants  in  equation  (1)  and 
the  RMS  deviation.  In  no  case  was  the  RMS  deviation  greater 
than  0.25  degrees  centigrade  and  in  most  cases  it  was  less 
than  0.05  degrees  centigrade. 

The  recording  instrument  used  was  an  Esterline  Angus  12- 
channel  Wheatstone  bridge  servo  recorder.  Each  channel  was 
sampled  every  fifteen  seconds.  A  chart  speed  of  8  inches  per 
hour  was  used.  The  charts  were  digitized  and  converted  to  tem¬ 
perature  before  analysis. 

The  recorder  was  calibrated  according  to  the  linear  equa¬ 
tion 

R  =  C  +  C  U  (2) 

o  1 

with  R  the  resistance  in  ohms,  C0  and  Cx  constants  and  U  the 
reading  in  recorder  units.  Resistance  between  about  2000  and 
6000  ohms  could  be  recorded  without  going  off-scale.  The  recor¬ 
der  was  calibrated  a  number  of  different  times  using  a  decade 
resistance  box.  A  least-square  method  was  used  to  determine 
the  constants  in  equation  (2)  for  each  calibration.  The  RMS 
values  ranged  between  2.5  and  7.0  ohms  for  the  various  calibra¬ 
tions.  Because  of  the  non-linear  relationship  between  tem¬ 
perature  and  resistance  in  the  probe,  a  small  uncertainty  in 
the  resistance  will  result  in  a  larger  uncertainty  in  tempera¬ 
ture  for  temperatures  near  25 °C  than  for  temperatures  near  0. 

An  RMS  deviation  of  7  ohms  in  resistance  corresponds  to  an  RMS 
deviation  of  0.09  degrees  for  temperatures  near  25°C  and  about 
0.027  degrees  for  temperatures  near  0. 
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Calibrations  performed  at  various  times  yielded  consistent 
values  for  the  constant  C ,  but  significantly  different  values 
for  the  constant  Cq .  This  implies  that  while  the  slope  of  re¬ 
sistance  versus  recorder  reading  remained  stable,  the  instru¬ 
mental  zero  shifted.  Since  at  least  one  of  the  twelve  channels 
was  always  shorted  out  with  a  constant  resistor  of  known  value, 
the  value  of  CQ  in  equation  (2)  could  be  determined  for  any  time 
between  calibrations. 

As  a  further  check  on  the  reliability  of  the  temperature 
data,  comparisons  were  made  with  23  BT's  taken  in  1969  by  the 
personnel  of  the  Naval  Underwater  Systems  Center,  Systems  Meas¬ 
urement  Branch  at  Dresden,  N.Y.  The  RMS  deviations  between 
temperatures  obtained  from  thermistors  and  from  BT's  ranged 
between  0.2  and  0.8°C  for  the  five  thermistors  in  use  at  the 
time.  When  a  second  set  of  five  thermistors  was  emplaced  190 
feet  from  the  original  string  in  1971,  discrepancies  were  noted 
between  pairs  of  thermistors  at  the  same  depths  (Figures  23 
through  27) .  It  was  determined  that  the  absolute  temperatures 
obtained  from  the  newer  thermistors  were  in  error  although 
relative  temperatures  from  any  single  thermistor  were  reliable. 

It  was  decided  that  the  difficulties  of  obtaining  a  good 
absolute  calibration  for  data  obtained  in  the  field  from  ther¬ 
mistors  which  had  shown  consistently  reliable  calibrations  in 
the  laboratory  were  due  to  leaks  in  both  the  material  of  the 
potting  compound  and  the  conducting  cable.  A  new  set  of  ther¬ 
mistors  was  built,  using  a  different  potting  compound  and 
heavier  cable.  These  thermistors  were  emplaced  in  September 
of  1971.  Comparison  with  BT's  showed  these  new  thermistors  to 
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be  much  more  reliable.  It  can  be  seen  from  Figures  31,  32, 

34,  36,  44,  46,  51  and  53  that  pairs  of  thermistors  at  simi¬ 
lar  depths  give  consistent  results. 

Temperature  profiles  at  both  barge  sites  were  obtained 
at  various  times  during  the  four  years  of  the  experiment  using 
a  bathythermograph.  The  instrument  used  had  a  temperature 
range  of  -2  to  +30 °C  and  could  record  to  a  depth  of  275  meters. 

A  temperature  versus  depth  profile  was  etched  on  a  glass  slide 
in  the  instrument  as  it  was  lowered  and  raised  through  the 
water.  The  glass  slide  was  then  removed  from  the  instrument, 
inserted  in  a  viewer  and  the  profile  was  redrawn  on  graph  paper. 
Experimental  set-up 

The  Naval  Underwater  Systems  Center,  Systems  Measurement 
Branch  (NUSC,  SMB)  at  Dresden,  has,  as  its  primary  function, 
the  testing  and  calibration  of  transducers  developed  and  used 
by  the  U.S.  Navy.  Seneca  Lake  was  chosen  as  the  site  for  this 
facility  because  of  its  great  depth.  Until  1970,  tests  were 
conducted  on  a  facility  composed  of  two  barges,  the  Transducer 
Calibration  Platform  (TCP)  and  the  WFNX-22  lighter,  which  were 
tied  together  and  moored  near  the  center  of  the  lake  at  site  1 
on  the  map  in  Figure  1. 

A  thermistor  string  with  five  probes  at  depths  of  25,  50, 
75,  100  and  150  feet  was  installed  in  August,  1968.  A  steel 
cable  with  a  60-pound  weight  at  the  end  was  used  to  hold  the 
string  vertical.  Each  probe  with  its  conducting  cable  was 
taped  to  the  steel  wire,  commencing  at  the  proper  depth  and 
continuing  up  to  the  surface.  The  support  wire  was  tied  to  a 
cleat  on  the  east  side  of  the  TCP  and  the  five  conducting  cables 
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Figure 


77° 


Location  map  of  Seneca  Lake  and  Facilities.  No. 
the  site  of  the  TCP  and  no.  2  is  the  site  of  the 


1  is 
SMP. 


The  scale  is  in  miles. 
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were  run  into  the  laboratory  on  the  barge  where  they  were  con¬ 
nected  to  the  recording  instrument.  Temperatures  were  conti¬ 
nuously  recorded  in  this  manner  until  the  beginning  of  February, 
1969,  when  the  instruments  were  removed.  The  thermistor  chain 
was  installed  again  in  May,  1969  and  removed  in  November  of 
that  year.  It  was  installed  once  more  in  May,  1970. 

In  August,  1970,  a  new  facility,  the  Systems  Measurement 
Platform  (SMP) ,  was  completed  and  the  NUSC,SMB  moved  their  base 
of  operations  to  it.  It  is  a  large,  floating  facility  composed 
of  two  barges  connected  by  a  superstructure  and  shaped  like  the 
letter  "H" .  It  is  permanently  moored  at  its  four  corners  at  a 
location  near  the  middle  of  the  lake  about  2.4  miles  north  of 
the  TCP  (see  map.  Fig.  1).  The  SMP  is  oriented  with  the  two 
hulls  parallel  to  each  other  and  the  bows  pointing  north.  The 
thermistor  string  was  hung  at  the  north  end  of  the  east  hull 
and  cables  were  run  along  the  deck  and  into  the  laboratory 
which  was  located  on  the  superstructure  connecting  the  barges. 

In  September,  1970,  a  second  thermistor  chain  was  installed  at 
the  south  end  of  the  east  hull  which  is  190  feet  long.  Origi¬ 
nally  consisting  of  three  probes  at  25,  100  and  150  feet  it  was 
completed  by  the  addition  of  two  more  probes  at  50  and  75  feet 
in  November.  Both  of  these  strings  were  removed  in  January, 
1971,  and  installed  again  in  May. 

Because  of  the  problems  of  obtaining  absolute  temperatures 
from  the  data  obtained  from  these  thermistors  (especially  those 
on  -the  south  chain)  a  new  set  was  built,  using  a  different  pot¬ 
ting  compound  and  heavier  cable  that  included  a  strength  member, 
eliminating  the  necessity  of  an  extra  wire.  Each  thermistor 
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was  hung  separately  with  its  own  15-pound  lead  weight.  This 
much  simplified  the  problem  of  checking  on  an  individual  ther¬ 
mistor  and  permitted  the  probes  to  be  lowered  to  the  proper 
depth  much  more  accurately  (on  the  old  strings,  the  steel  wire 
with  the  weight  had  a  tendency  to  slip  past  the  thermistor 
cables  as  the  string  was  lowered  into  the  water) .  These  new 
probes  were  installed  at  the  two  ends  of  the  east  hull  of  the 
SMP  at  the  same  depths  as  the  old  ones,  in  September,  1971.  In 
October,  a  sixth  probe,  consisting  of  an  old  thermistor  con¬ 
nected  to  the  new,  heavier  cable,  was  hung  from  the  north  end 
of  the  SMP  at  a  depth  of  225  feet.  The  installation  was  again 
disassembled  and  removed  in  December. 

The  instruments  were  removed  every  winter  because  the  lake 
is  nearly  isothermal  throughout  most  of  the  winter  and  early 
spring,  so  very  little  useful  data  would  be  obtained  with  the 
instruments  used.  During  this  time  the  thermistors  and  the 
recorder  were  calibrated  and  if  it  needed  to  be,  the  recorder 
was  serviced. 

In  the  summer  of  1971,  thermistors  were  installed  at  site 
1  (Figure  1)  with  the  purpose  of  measuring  the  speed  of  propa¬ 
gation  of  internal  waves  between  the  two  sites.  In  July  a 
Rustrak  2-channel  temperature  recorder  was  installed  on  the 
TCP.  The  thermistor  depths  were  50  and  100  feet.  In  laboratory 
calibrations  it  had  become  evident  that  the  thermistors  used 
with  this  recorder  had  considerable  drift.  The  data  obtained 
were  therefore  used  only  to  obtain  times  for  internal  waves. 

This  recorder  was  in  operation  for  three  weeks  in  July  and 
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August,  after  which  the  TCP  was  removed  from  this  site  and  the 
lighter  was  moored  there.  For  one  week  in  the  middle  of  Sep¬ 
tember  and  then  from  late  October  until  December,  temperatures 
were  recorded  at  site  1  using  a  thermistor  arrangement  similar 
to  that  employed  on  the  SMP.  The  recording  instrument  was  an 
Esterline-Angus  24-point  chart  recorder  with  a  scale  running 
from  0  to  1  millivolt.  The  thermistors  were  connected  through 
bridge  circuits  to  provide  voltage  outputs. 

Because  the  thermistors  used  were  old  ones  and  because  of 
problems  with  the  bridge  circuits  used  (there  was  some  crossover 
between  channels)  the  data  cannot  be  used  to  obtain  reliable 
temperatures.  However,  internal  waves  can  be  seen  and  used  to 
determine  propagation  velocities  between  the  two  sites. 

In  Figure  2  are  shown  all  the  intervals  when  useful  data 
was  obtained  during  the  four  years  the  instruments  were  in 
operation . 

Data  handling 

In  1968  data  was  digitized  from  mid-August  to  the  end  of 
November.  Each  record  had  at  least  two  time  marks  that  had  been 
put  on  by  hand  while  the  data  was  being  recorded.  From  the  time 
marks  and  the  known  chart  speed,  marks  were  put  on  the  record 
at  hourly  intervals.  Values  were  then  read  off  the  record  for 
each  of  the  five  probes  and  for  one  channel  shorted  with  a  re¬ 
sistor.  These  were  then  punched  on  IBM  cards  and  eventually 
put  on  tape.  The  data  is  continuous  except  for  a  few  gaps  when 
the  paper  on  the  recorder  wasn't  changed  in  time  and  for  some 
power  failures.  The  longest  segment  of  no  data  was  in  mid- 
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November  and  amounted  to  69  hours.  Eight  segments  of  data  at 
times  of  short-period  internal  wave  activity  were  later  digi¬ 
tized  directly  from  the  records  to  IBM  cards  using  a  digitizing 
machine.  These  segments  were  digitized  at  one-minute  intervals. 

The  1969  data  from  May  to  November  was  digitized  directly 
from  the  records  using  the  digitizing  machine.  Rather  than  di¬ 
gitizing  at  any  fixed  time  interval,  a  point  was  digitized  on 
a  channel  whenever  a  change  in  slope  occurred.  The  times  of 
all  the  digitized  points  was  later  calculated  by  computer  for 
all  records  that  contained  at  least  two  time  marks.  Unfortu¬ 
nately,  there  are  a  few  records  from  1969  that  contain  only  one 
or  no  time  marks.  Thus,  while  the  records  digitized  span  an 
interval  of  188  days  between  May  8  and  November  5,  there  are 
less  than  108  days  of  usable  data.  Of  the  81  days  of  missing 
data,  56  occur  in  May,  June  and  July. 

The  data  obtained  during  1970  was  similar  to  that  obtained 
in  the  previous  two  years  and  was  not  digitized.  The  1971  data 
was  digitized  only  during  times  of  short-period  internal  wave 
activity  using  the  same  method  employed  with  the  1969  data.  Be¬ 
cause  of  difficulties  encountered  with  the  recorder,  there  were 
only  a  few  short  intervals  of  reliable  data  between  July  30  and 
October  20. 

Once  digitized,  the  data  was  converted  into  resistance 
using  equation  2  and  the  reading  from  a  channel  shorted  by  a 
resistor.  Equation  1  was  used  to  convert  the  data  to  temperature. 
The  final  form  of  the  data  are  time  series  of  temperature  for 
specific  depths. 
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GENERAL  FEATURES  OF  THE  TEMPERATURE  DATA 
Yearly  temperature  cycle 

Temperature  at  all  levels  in  the  lake  fluctuates  over  a 
large  range  of  time  scales,  from  climatic  changes  extending 
over  centuries  to  turbulence  with  periods  of  only  a  few  seconds. 
Two  prominent  time  scales  in  Seneca  Lake  are  the  seasonal  one 
and  one  due  to  internal  surges  with  periods  of  ten  minutes  to 
a  few  days . 

Seasonal  changes  in  the  temperature  structure  of  lakes 
have  been  observed  and  discussed  by  many  authors.  A  comprehen¬ 
sive  summary  and  bibliography  was  given  by  Hutchison  (1957) . 

In  his  classification,  Seneca  Lake  should  be  temperate  or  dimic- 
tic,  circulating  twice  a  year.  In  late  autumn  or  early  winter 
and  again  in  mid-spring,  the  lake  is  isothermal  at  the  tempera¬ 
ture  of  maximum  density,  4°C  (39°F) ,  and  mixing  throughout.  In 
between,  the  surface  water  should  cool  (and  freeze  if  the  winter 
were  long  enough  or  cold  enough) ,  the  lake  should  become  strati¬ 
fied  with  a  negative  temperature  gradient  and  should  not  mix 
until  the  surface  warmed  to  4°C.  However,  the  few  winter  pro¬ 
files  obtained,  although  showing  a  weak  negative  temperature 
gradient,  indicate  that  Seneca  Lake  does  not  conform  to  this 
simple  model.  Figure  3  shows  the  typical  temperature  structure 
of  the  lake  at  various  times  of  the  year.  The  curves  shown 
were  selected  from  the  bathythermograph  records. 

In  late  autumn  or  early  winter  the  lake  is  isothermal  above 
4°C  (Figure  3,  December  31,  1969).  Later  in  the  winter  there 
is  an  inverse  temperature  structure  (Figure  3,  February  16,  1968) 


DEPTH  r.  r-_~  DEPTH  IN  FEET  DEPTH  IN  PEET 


15 


LTD-  31.  19G9 


FEB-  1G.  19G8 


TEMP  IN  DEG  CENT 
MAY  B>  1970 


.JUNE  25.  1968 


A°R-  17,  1969 

T 


TEMP  IN  DEb  CENT 
AUG-  15,1968 


SEPT-  25,  1969 


NOV-  10,  1971 


MDV-  20,  1969 


Figure  3.  Vertical  temperature  structure  of  Seneca  Lake  at  various 
times  of  the  year.  Bathythermographs  from  1968,  1969 
and  1970  were  taken  at  barge  site  1;  from  1971  at  barge 


site  2 
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but  the  surface  layer  is  very  thick  (over  200  feet)  and  not 
much  colder  than  the  deep  layer.  By  early  spring,  the  lake  is 
again  isothermal  but  below  4°C  (Figure  3,  April  17,  1969). 

This  pattern  of  the  lake  becoming  isothermal  at  a  temperature 
above  4°C  in  early  winter,  cooling  and  developing  two  thick 
layers  with  a  very  small  temperature  difference  in  mid-winter 
and  becoming  isothermal  below  4°C  in  late  winter  and  early 
spring  is  seen  in  other  winter  profiles.  Berg  (1963)  in  dis¬ 
cussing  nearby  Cayuga  Lake,  which  is  very  similar  to  Seneca 
Lake,  maintains  that  except  during  unusually  cold  winters  that 
lake  circulates  completely  throughout  the  winter.  He  refers  to 
data  obtained  by  other  workers,  especially  Henson  et  al.  (1961), 
to  support  his  conclusion  and  proposes  that  the  term  warm  mono- 
mi  ctic,  with  a  change  in  the  definition  given  by  Hutchinson 
(1957)  ,  be  used  to  classify  lakes  such  as  Cayuga  (and  possibly 
Seneca)  that  have  one  circulation  period  annually. 

It  is  uncertain  whether  Seneca  Lake  should  be  considered 
as  dimictic  or  warm  monomictic  in  Berg's  sense.  The  inverse 
temperature  structure  seen  in  the  few  profiles  taken  in  mid¬ 
winter  indicate  that  the  lake  is  dimictic  but  the  stratification 
is  weak.  In  any  case,  the  lake  is  isothermal  and  freely  circu¬ 
lating  for  a  good  part  of  the  winter  and  early  spring  with  tem¬ 
peratures'  both  above  and  below  4°C.  Seneca  Lake  may  be  a  bor¬ 
derline  case  between  dimictic  and  monomictic;  some  winters  the 
inverse  stratification  may  be  strong  and  the  lake  would  behave 
more  like  the  classical  dimictic  model,  while  in  other  years 
the  stratification  might  not  develop  at  all  and  the  lake  would 
be  monomictic. 
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In  late  April  or  May,  the  surface  warms,  the  lake  becomes 
stratified  and  a  thermocline  begins  to  develop  (Figure  3,  May 
8,  1970).  By  the  end  of  June  the  lake  consists  of  two  isothermal 
layers,  a  warm  upper  layer,  the  epilimnion,  and  a  cold  lower 
layer,  the  hypolimnion,  separated  by  a  strong  thermocline  with 
an  approximately  exponential  temperature  gradient  (Figure  3, 

June  25,  1968).  Throughout  most  of  the  summer,  the  epilimnion 
warms,  reaching  a  maximum  temperature  in  August  (Figure  3,  Aug¬ 
ust  15,  1968)  and  then  begins  to  cool  (Figure  3,  September  25, 
1969) .  Throughout  most  of  the  fall  the  upper  layer  continues 
to  cool  and  thicken  (Figure  3,  November  10,  1971  and  November 
20,  1969).  By  late  December  or  early  January,  the  lake  is  once 
again  isothermal  and  mixing  throughout.  The  lake  need  not  be 
at  4°C  when  this  occurs.  Berg  (1963)  points  out  that  Cayuga 
Lake  becomes  isothermal  at  a  temperature  appreciably  above  4°C. 
Heat  budget 

Another  way  of  looking  at  the  yearly  temperature  cycle  in 
a  lake  is  to  consider  the  heat  budget.  A  mean  temperature  can 
be  determined  for  the  lake  at  any  time  that  the  thermal  struc¬ 
ture  is  known.  The  mean  temperature  is  the  temperature  that  the 
lake  would  be  at  if  it  were  to  be  completely  mixed  without  any 
gain  or  loss  of  heat.  To  determine  the  mean  temperature  the 
lake  is  divided  into  a  number  of  layers,  the  average  temperature 
in  each  layer  is  calculated  and  multiplied  by  the  volume  of 
that  layer.  The  results  for  all  the  layers  are  summed  up  and 
divided  by  the  total  lake  volume.  Birge  and  Juday  (1914)  made 
the  first  attempt  to  obtain  mean  temperatures  and  heat  budgets 
in  the  Finger  Lakes  from  data  collected  from  1910  through  1912. 
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They  obtained  temperatures  in  Seneca  Lake  on  five  days  at 
the  beginning  of  August,  1910  and  on  February  10  and  September 
1,  1911.  Data  was  also  obtained  on  September  5,  1914  and  on 
August  29,  1918.  Only  a  single  temperature  profile  was  used 
for  each  determination  of  a  mean  temperature.  A  mean  tempera¬ 
ture  determined  in  this  fashion  will  be  valid  only  if  the  tem¬ 
perature  profile  used  is  representative  of  the  average  condi¬ 
tions  of  the  whole  lake  and  not  only  of  the  point  where  it  was 
obtained.  For  this  to  be  the  case,  either  the  isotherms  in 
the  lake  cannot  be  tilted  or  if  they  are,  the  temperature  pro¬ 
file  must  be  taken  at  a  point  in  the  lake  where  "average"  con¬ 
ditions  prevail.  Birge  and  Juday  were  well  aware  of  internal 
waves  in  lakes  (they  called  them  "temperature  seiches")  and 
endeavored  to  make  their  measurements  near  the  "center  of  oscil¬ 
lation"  of  the  lake.  They  obtained  their  measurements  near  the 
center  of  Seneca  Lake,  assuming  that  the  internal  waves  were 
linear,  first  mode  oscillations  with  maximum  amplitude  at  the 
ends  of  the  lake  and  a  node  near  the  center.  It  will  become 
evident  from  the  data  presented  in  this  paper  that  the  internal 
waves  in  Seneca  Lake  are  non-linear  and  travel  through  the  lake 
as  a  surge.  Therefore,  temperature  profiles  obtained  even  near 
the  center  of  the  lake  may  not  represent  a  true  picture  of  the 
complete  thermal  structure. 

Figure  4  shows  the  variations  in  the  mean  temperature  of 
the  lake  over  the  course  of  three  different  years.  The  data 
is  from  BT's  taken  at  site  1  in  1968  and  1969  and  site  2  in 
1971  (Figure  1) .  Again  it  must  be  emphasized  that  because  of 
the  distortions  of  the  thermal  structure  caused  by  internal 
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TIME  IN  DAYS  EROM  JAN*  1  OOOO  EST 


Figure  4.  Mean  temperature  of  Seneca  Lake  for  three  different 

years.  The  vertical  dashed  lines  separate  the  months. 
The  solid  curve  is  for  1968.  The  dashed  curve  with  a 
long  line  and  a  short  space  is  for  1969.  The  dashed 
curve  with  a  long  and  short  line  separated  by  a  short 
space  is  for  1971. 
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surges,  any  individual  mean  temperature  plotted  in  Figure  4 
must  not  be  considered  very  reliable  (for  example,  the  tempera¬ 
tures  plotted  for  July  14  and  30,  1971;  the  former  is  probably 
too  high  while  the  latter  is  too  low).  However,  the  general 
pattern  of  the  lake's  thermal  cycle  can  be  seen  in  the  figure. 
The  lake  warms  throughout  the  spring  and  early  summer,  the  tem¬ 
perature  increasing  most  rapidly  in  June.  By  the  end  of  July 
the  temperature  is  increasing  at  a  very  moderate  rate,  reaching 
a  peak  in  September  or  October.  There  is  then  a  moderate  dec¬ 
line  in  temperature  followed  by  a  more  rapid  loss  of  heat  in 
November  and  December.  This  general  pattern  is  seen  for  all 
three  years  plotted.  The  1968  and  1969  data  are  very  similar 
but  the  1971  data  appears  erratic  throughout  the  summer  and 
early  fall.  For  comparison,  Birge  and  Juday  (1922)  obtained 
summer  temperatures  ranging  from  7.34  to  8.27°C  from  data  ob¬ 
tained  during  four  different  years. 

The  data  shown  in  Figure  4  can  be  represented  in  a  few 
other  ways.  If' the  mean  temperature  is  multiplied  by  the  total 
volume  of  water  in  the  lake,  the  total  heat  content  in  calories 
above  0°C  will  be  obtained.  This  is  a  less  useful  figure  than 
the  mean  temperature  because  it  cannot  readily  be  compared  with 
other  lakes  (the  heat  content  of  a  large,  cold  lake  can  be 
greater  than  that  of  a  smaller,  warmer  one) .  However,  if  we 
then  divide  the  total  heat  content  by  the  surface  area  of  the 
lake,  we  obtain  the  Birgean  heat  budget,  ©b,  which  is  the  heat 
per  unit  surface  area  necessary  to  bring  the  lake  from  0°C  to 
its  mean  temperature.  Since  almost  all  of  the  heat  gains  and 
losses  of  the  lake  occur  at  its  surface,  this  is  a  very  useful 


21 


figure.  Subtracting  the  Birgean  heat  budget  at  the  time  of 

the  lowest  mean  temperature  in  winter  from  that  of  the  highest 

in  summer,  results  in  the  annual  heat  budget, 6  ,  which  is 

.D  a 

the  total  amount  of  heat  per  unit  area  that  must  enter  the 
lake  as  it  warms  or  leave  as  it  cools.  The  summer  heat  income, 

©ks  ,  is  the  amount  of  heat  per  unit  area  that  enters  the  lake 
between  the  time  when  the  lake  is  isothermal  at  4°C  and  when 
it  is  at  its  warmest.  The  term  wind-distributed  heat  is  some¬ 
times  used  for  this  quantity.  Because  the  lake  is  normally 
stratified  above  4°C  and  solar  radiation  can  only  heat  the 
upper  few  meters  of  water  directly,  most  of  this  heat  is  deli¬ 
vered  to  the  lower  layers  of  the  lake  by  mixing  of  water  with 
the  energy  coming  from  the  wind. 

Birge  and  Juday  (1914)  obtained  only  one  winter  mean  tem¬ 
perature  for  the  lake,  3.39°C  on  February  10,  1911.  From  the 
summer  temperatures  measured  in  1910  and  1911,  they  calculated 

2 

two  annual  heat  budgets  and  obtained  38,300  and  35,100  cal. /cm. 
Two  winter  mean  temperatures  have  been  calculated  from  BT  data; 
3.52°C  on  February  16,  1968  and  3.34°C  on  April  17,  1969.  As¬ 
suming  the  maximum  summer  temperature  to  be  between  8.0  and 
8.5°C  and  the  minimum  winter  temperature  between  3.2  and  3.5°C, 
the  annual  heat  budget  will  be  between  40,000  and  47,000  cal./ 
cm.  for  most  years  (the  results  of  Birge  and  Juday  are  low 
because  the  summer  temperatures  they  used,  7.71°C  in  1910  and 
7.35°C  in  1911  are  low;  they  probably  are  not  the  maximum  tem¬ 
peratures  achieved  by  the  lake  during  those  summers) .  The 
summer  heat  income  for  Seneca  Lake  was  estimated  to  be  between 
about  35,500  and  almost  40,000  cal. /cm. ^ 


The  four  values 
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obtained  by  Birge  and  Juday  (1922)  ranged  from  less  than 
30,000  to  almost  39,000  cal. /cm. ^ 

A  simple  calculation  shows,  that  if  we  assume  that  the 
maximum  and  minimum  temperatures  occur  six  months  apart,  the 
average  heat  flow  through  the  lake  surface  is  between  220  and 
255  cal. /cm.  per  day.  From  Figure  4  it  is  seen  that  the  heat 
flow  through  the  surface  is  not  constant  but  varies  through¬ 
out  the  year.  It  is  maximum  in  the  spring  and  fall  and  near 
zero  in  summer  and,  presumably,  in  winter.  The  heat  flow 
during  any  time  interval  between  two  observed  mean  temperatures 
can  also  be  calculated.  Unfortunately,  because  of  the  large 
error  bounds  of  the  individual  temperature  points,  a  graph  of 
heat  flow  versus  time  would  look  very  erratic.  Instead,  the 
curves  in  Figure  4  were  smoothed  by  eye  and  a  composite  smooth 
curve,  representing  an  "average"  year  was  drawn.  This  curve 
was  then  used  to  determine  heat  flow  for  various  months.  The 
heat  gain  is  greatest  during  June,  exceeding  400  cal. /cm. 2  per 
day.  This  can  be  compared  with  the  total  solar  radiation  in 
June  at  the  latitude  of  the  lake  which  is  about  600  cal. /cm. 2 
per  day.  This  value  is  augmented  by  the  scattered  light  of  the 
sky  and  is  diminished  whenever  the  sun  is  obscured  by  clouds. 

The  heat  flow  diminishes  throughout  the  summer  reaching 
zero  in  September,  after  which  there  is  a  heat  loss  that  con¬ 
tinues  to  grow  throughout  the  fall.  It  is  over  300  cal. /cm. 2 
per  day  in  November.  Because  of  the  paucity  of  winter  data, 
the  time  and  value  of  the  maximum  heat  loss  could  not  be  esti¬ 
mated. 

In  addition  to  the  total  heat  in  the  lake,  the  heat  in 


23 


each  ten-meter  layer  has  been  calculated  for  a  number  of  points 
in  1968,  1969  and  1971.  Figure  5  shows  the  total  heat  in  the 
upper  20  meters  of  the  lake  for  the  three  years.  It  can  be 
seen  that,  like  the  total  heat  (Figure  4) ,  the  heat  content  of 
the  upper  layer  increases  rapidly  for  most  of  the  spring  and 
early  summer,  reaching  a  maximum  in  August  or  September.  The 
BT  records  (Figure  3)  show  the  same  behavior;  that  most  of  the 
heat  entering  the  lake  during  this  period  remains  in  the  upper 
20  meters.  In  October  and  early  November,  however,  at  a  time 
when  the  lake  as  a  whole  is  in  thermal  equilibrium  or  losing 
heat  at  a  slow  rate,  the  upper  layer  is  rapidly  losing  heat. 
Figure  6,  in  which  the  total  heat  between  20  and  40  meters  has 
been  plotted,  shows  that  some  of  the  heat  lost  by  the  upper  20 
meters  has  been  transported  down.  This  layer  gains  heat  slowly 
throughout  the  spring  and  early  summer  until  September  or  Octo¬ 
ber  when  there  is  a  rapid  heat  gain  followed  by  a  rapid  heat 
loss.  Figure  7  shows  the  total  heat  in  the  40  to  70  meter 
layer.  In  this  layer  there  is  almost  no  heat  gain  from  June 
through  October.  In  October,  the  layer  begins  to  warm,  the 
rate  becoming  maximum  in  November,  after  which  the  layer  cools. 
A  similar  pattern  is  seen  in  the  bottom  region  of  the  lake. 
Figure  8  shows  the  total  heat  between  70  meters  and  the  bottom 
of  the  lake  at  188.4  meters. 

In  summary,  the  five  sets  of  curves  plotted  in  Figures  4 
through  8  and  the  yearly  BT  cycle  (Figure  3)  illustrate  the 
annual  thermal  cycle  of  the  lake.  In  May  and  early  June,  the 
upper  region  of  the  lake  warms,  with  an  isothermal  upper  layer, 
the  epilimnion,  and  a  thermocline  forming.  The  epilimnion  is 
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Figure  5.  Total  heat  above  0°C  between  the  surface  and  20  meters 
in  Seneca  Lake  during  the  years  1968,  1969  and  1971. 


i 
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Figure  6. 


Total  heat  above  0°C  between  20  and  40  meters  in  Seneca 
Lake  during  the  years  1968,  1969  and  1971. 
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Figure  7. 


Total  heat  above  0°C  between  40  and  70  meters  in  Seneca 
Lake  during  the  years  1968,  1969  and  1971. 


27 


Figure  8.  Total  heat  above  0°C  between  70  meters  and  the  deepest 
part  of  Seneca  Lake  (188.4  meters)  during  the  years 


1968,  1969  and  1971. 
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confined  within  the  upper  20  meters  while  the  water  below  70 
meters  feels  little  of  the  effects  of  this  warming  of  the  lake. 
When  the  epilimnion  is  at  its  highest  temperature,  most  of  the 
heat  that  has  been  added  to  the  lake  is  still  confined  to  the 
upper  70  meters.  It  is  only  after  the  epilimnion  begins  to 
cool  that  it  thickens,  and  considerable  amounts  of  heat  are 
carried  to  the  lower  regions  of  the  lake. 
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LONG  PERIOD  INTERNAL  WAVES 

Introduction 

An  examination  of  the  data  reveals  that,  in  addition  to 
the  yearly  temperature  cycle,  there  are  large-scale  temperature 
oscillations  with  periods  of  tens  of  hours  to  a  few  days.  Fi¬ 
gure  9  shows  the  temperature  at  five  levels  in  the  lake  from 
mid- August  to  the  end  of  November,  1968.  The  data  was  digitized 
at  one-hour  intervals.  In  many  instances  it  is  seen  that  a 
minimum  temperature  recorded  at  one  depth  is  below  the  preceding 
and  succeding  maximum  temperatures  recorded  by  the  next  deepest 
thermistor.  Since  the  top  four  thermistors  were  spaced  at  25- 
foot  intervals,  this  implies  water  motions  with  vertical  ampli- 
tud  of  at  least  12-1/2  feet.  There  are  a  few  cases  where  a 
maximum  or  minimum  temperature  moves  past  what  would  be  consi¬ 
dered  the  mean  temperature  of  the  thermistor  above  or  below  it, 
indicating  25- foot  amplitudes. 

The  oscillations  seen  in  Figure  9  cannot  be  attributed  to 
actual  changes  in  the  thermal  structure  of  the  lake  as  a  whole 
because  the  periods  seen  are  much  too  short.  It  was  shown  in 
the  sections  on  the  seasonal  changes  in  the  lake  and  on  the 
heat  budget  that  these  changes  have  time  scales  of  weeks  and 
months  and  not  days.  It  is  therefore  evident  that  these  large 
temperature  fluctuations  recorded  by  the  thermistors  cannot  be 
due  to  warming  and  cooling  of  water  but  must  be  due  to  different 
water  masses  of  fairly  constant  temperature  moving  past  the 
thermistors,  which  is  one  of  the  criteria  of  wave  motion.  The 
motion  also  appears  to  be  periodic,  although  the  period  is  not 
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constant.  Other  criteria  of  wave  motion  such  as  energy  trans¬ 
port  cannot  be  determined  from  measurements  made  at  one  location. 
Because  they  have  been  seen  in  other  lakes,  it  will  be  assumed 
here  that  the  temperature  oscillations  are  due  to  internal  waves 
although  some  of  the  motion  is  undoubtedly  turbulent. 

Internal  waves  in  narrow  lakes  have  been  well  known  for 
more  than  half  a  century.  Most  of  the  early  investigators 
(Watson,  1904;  Wedderburn,  1907,  1909,  1911  and  1912;  Wedder- 
burn  and  Williams,  1911)  could  not  obtain  continuous  recordings 
of  temperature  and  thus  could  only  analyze  internal  waves  with 
periods  of  hours  or  days.  Mortimer  (1953)  gives  a.  good  summary 
of  the  work  done  up  to  the  middle  of  this  century,  showing  that 
the  assumption  of  internal  seiches  in  two-layered  basins  gives 
good  agreement  between  theory  and  the  observed  data  in  most 
cases.  (Henson  (1959)  gives  evidence  of  internal  seiches  with 
a  57-hour  period  in  nearby  Cayuga  Lake  in  August,  1951.  The 
amplitude  was  over  20  meters  at  his  station  which  was  about 
10  miles  from  the  southern  end  of  the  lake  which  is  38  miles 
long.)  It  would  appear  that  the  first  attempt  to  understand 
the  long-period  oscillations  observed  in  Seneca  Lake  should  be 
along  these  lines . 

Theory  of  internal  waves  in  a  two-layered  fluid 

Internal  waves  can  be  defined  as  waves  in  which  the  maxi¬ 
mum  amplitude  occurs  within  the  fluid  rather  than  at  the  sur¬ 
face.  They  are  an  inherent  characteristic  of  a  stratified 
fluid  that  has  been  disturbed  from  equilibrium  and  are  analogous 
to  surface  waves  (indeed,  surface  waves  can  be  considered  to  be 
a  special  case  of  a  two-layered  fluid  with  the  density  of  the 
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upper  layer  taken  as  equal  to  zero) .  A  period  equation  for 
internal  waves  in  a  two-layered  fluid  can  be  derived  (see,  for 
example,  Neumann  and  Pierson,  1966).  In  the  derivation,  the 
Coriolis  force  has  been  ignored;  this  is  justified  for  narrow 
bodies  of  water  like  Seneca  Lake. 

w4(p  coth (kh) coth (kh7 ) +p 7 ) -  w2kpg ( coth (kh) +coth (kh 7 ) ) 

+  g2k2 ( p-p 7 )  =  0  ( 3 ) 

In  equation  3,  w  is  the  frequency,  k  the  wave  number  (  =  2tt/A, 
with  A  the  wave  length) ,  p  the  density,  h  the  thickness  of  the 
layer  and  g  the  acceleration  of  gravity.  The  primed  symbols 
refer  to  the  upper  layer,  the  unprimed  to  the  lower. 

Equation  3  is  quadratic  in  w 2  and  can  easily  be  solved 
for  w2  in  terms  of  k.  However,  for  waves  of  wavelengths  on 
the  order  of  the  length  of  Seneca  Lake,  the  expressions  kh  and 
kh 7  will  be  very  close  to  zero  and  the  terms  coth (kh)  and 
coth(kh7)  can  be  replaced  by  (kh)-^  and  (kh7)”^-.  Making  this 
simplification  and  replacing  w  by  kc  (c  is  the  phase  velocity) , 
equation  3  can  be  rewritten: 

c4-c2g (h+h  7 )  +  g2hh  7  (P.~.P')  =  0  (4) 

P 

This  equation  has  two  roots  which  are,  approximately: 

Cl2  =  g (h+h ' ) 

c 2 2  “  g  “f-d  -_£!  )  (5) 

h+h  p 

The  first  root  is  the  familiar  long  wave  velocity  of  a  surface 
gravity  wave.  The  second  root  corresponds  to  the  internal 
wave.  The  phase  velocity  of  the  surface  waves  on  Seneca  Lake 
(whose  mean  depth  is  290  feet)  will  be  about  100  feet  per  second. 
For  the  internal  waves  the  phase  velocity  will  be  on  the  order 
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of  1  foot  per  second  for  the  summer  and  early  fall  stratifi¬ 
cation,  the  exact  value  depending  on  the  lake  structure  at  the 
time . 

An  expression  for  the  ratio  of  the  amplitudes  of  an  inter¬ 
nal  wave  at  the  interface  and  the  surface  can  also  be  deter¬ 
mined: 

n  -  _  P (h+h ' )  (6) 

ri  '  h  ( p-p  1 ) 

with  n  the  amplitude  at  the  interface  and  n '  at  the  surface. 

The  conditions  in  Seneca  Lake  yield  values  of  this  ratio  on  the 
order  of  1000.  Thus,  internal  waves  with  amplitudes  of  tens  of 
feet  can  have  associated  surface  waves  with  amplitudes  of  less 
than  an  inch. 

The  natural  period  of  the  lake,  also  called  the  seiche 
period,  can  be  determined  for  both  the  surface  and  internal 
waves.  The  seiche  period  of  order  (or  mode)  n  is  that  period 
that  will  allow  exactly  n/2  wavelengths  to  fit  the  length  of 
the  lake: 

(-j-)A  =  L 

with  A  the  wavelength  and  L  the  length  of  the  lake.  The  period 
can  be  determined  from: 


T  = 


A 


2L  1 
n  c 


Substituting  for  the  phase  velocity,  c,  from  equation  5,  the 
surface  and  internal  wave  seiche  periods  are: 

2L 

| gtn+n- ;  \ 

(7) 


T  i = 


n 


T  = 

2 


[g(h+h'fj 
& 


2L  F  p  (h+h' )  •"] 
n  ]  ghh'  (p-p  ;  )J 
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The  fundamental  seiche  period  (n=l)  is  about  one  hour  in  Lake 
Seneca  for  the  surface  wave.  For  the  internal  wave  it  will  be 
between  45  and  95  hours,  depending  upon  the  thermal  structure 
at  the  time . 

The  non-linear  characteristics  of  the  observed  long-period 

waves 

It  is  at  this  point  that  the  earlier  investigators  would 
measure  the  period  of  the  internal  waves  and  compare  it  with 
the  period  predicted  from  equation  7.  Mortimer  (1953)  did  this 
for  a  number  of  lakes  and  showed  that  the  observed  period  was 
very  close  to  the  theoretical  uninodal  seiche  period  in  most 
cases  (see  his  Tables  1  and  2) .  This  approach  presents  pro¬ 
blems  when  it  is  attempted  with  the  data  obtained  from  Seneca 
Lake.  Figure  10  shows  the  temperature  data  obtained  during 
October  1968.  This  figure  is  similar  to  Figure  9  but  only  30 
instead  of  110  days  of  data  have  been  plotted  with  an  appropri¬ 
ate  expansion  of  the  time  axis. 

It  is  clearly  seen  that  there  are  no  sinusoidal  oscil¬ 
lations  with  periods  of  a  few  days.  There  are  a  number  of  large 
amplitude  variations  in  which  temperatures  decline  slowly  for 
some  tens  of  hours  and  then  return  to  their  former  level  within 
an  hour  (see  arrows  in  Figure  10) .  It  will  be  shown  later  in 
this  report  that  this  rapid  increase  in  temperature  is  accom¬ 
panied  by  large- amplitude  oscillations  with  periods  of  tens  of 
minutes  (which  do  not  appear  in  the  data  presented  in  Figures 
9  and  10  because  a  digitizing  interval  of  one  hour  has  been 
used.  They  can,  however,  be  seen  in  Figures  18  through  22) 
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EST.  The  arrows  mark  times  when  high  frequency  oscillations  were  re 
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and  that  they  are  non-linear  phenomena. 

Before  discussing  non-linear  waves  it  will  be  instructive 
to  look  at  temperature  data  that  has  been  filtered  to  reduce 
the  amplitude  of  short-period  oscillations.  Figure  11  shows 
the  1968  and  1969  data  after  having  been  put  through  cosine 
filters  with  widths  of  12  and  25  hours.  The  cosine  filter  re¬ 
duces  the  energy  of  those  oscillations  with  a  period  equal  to 
the  filter  width  by  50%.  Periods  twice  the  filter  width  are 
reduced  only  15%  while  periods  half  the  filter  width  are  com¬ 
pletely  lost.  It  is  seen  that  the  filtered  data  looks  smoother 
and  less  non-linear  than  the  original  data.  The  data  filtered 
with  a  12-hour  filter  (Figure  11,  bottom)  is  particularly  inter¬ 
esting.  If,  instead  of  being  recorded  every  15  seconds  (the 
sampling  interval  of  the  Esterline-Angus  recorder) ,  the  temper¬ 
ature  had  been  recorded  every  three  hours  (for  instance  by 
lowering  thermometers,  a  technique  used  before  thermistors  had 
been  developed)  and  a  smooth  curve  drawn  between  the  points,  the 
result  would  look  very  much  like  this  graph.  The  filtered  data 
also  show  that  the  upper  layers  of  the  lake  cooled  more  rapidly 
in  the  fall  of  1969  than  they  did  in  the  previous  fall.  This 
was  also  seen  in  Figure  5. 

Power  spectra  were  obtained  for  the  1969  hourly  data  and 
are  shown  for  the  five  thermistor  depths  in  Figures  12  through 
16.  There  are  very  few  peaks  that  depart  from  the  spectra  back¬ 
ground  by  more  than  the  90%  confidence  band  and  those  that  do, 
just  barely  do  so  and  in  many  cases  are  not  sharp  peaks.  It  is 
apparent  that  there  are  no  predominant  frequencies  in  the  data. 
This  was  also  seen  in  Figures  9  and  10  which  showed  large  tern- 
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Figure  12.  Power  spectrum  of  the  1969  hourly  data.  The  data  were 

from  May  8  to  November  12.  The  90%  confidence  band 
is  shown  in  the  upper  right  hand  corner.  Thermistor 
depth  is  25  feet. 
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Figure  13.  Similar  to  Figure  12.  Thermistor  depth  is  50  feet 
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DEPTH  =  75  FEET 


Figure  14. 


Similar  to  Figure  12.  Thermistor  depth  is  75  feet. 
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Similar  to  Figure  12. 


Thermistor  depth  is  100  feet. 
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Figure  16. 


Similar  to  Figure  12.  Thermistor  depth  is  150  feet. 
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perature  fluctuations  with  periodicities  of  a  few  days  but  the 
periods  were  not  constant  and  there  were  many  other  frequencies 
superimposed. 

The  most  striking  aspect  of  the  large-scale  temperature 
fluctuations  seen  in  Figure  10  is  its  asymmetry.  Mortimer 
(1955)  measured  temperatures  in  Loch  Ness  in  1952  by  means  of 
a  thermistor  chain.  He  found  that  the  main  internal  seiche 
was  characterized  by  a  sudden  rise  in  temperature  similar  to 
that  seen  in  Seneca  Lake.  Thorpe  (1971)  made  temperature  obser¬ 
vations  at  a  different  point  in  Loch  Ness  in  1970  and  confirmed 
the  asymmetrical  shape  of  the  seiche.  He  ruled  out  as  a  cause 
of  the  asymmetry  the  effects  of  the  thermocline  and  its  micro¬ 
structure,  the  breaking  of  the  seiche  wave  as  it  approaches 
the  end  of  the  Loch  and  the  effects  of  the  wind.  It  was  his 
conclusion  that  the  main  cause  of  the  asymmetry  observed  in 
Loch  Ness  was  non-linear  effects.  To  show  this,  he  introduced 
an  expansion  parameter,  6  (which  is  a  function  of  the  amplitude 
of  the  seiche  motion,  the  length  of  the  lake,  and  the  thickness 
and  densities  of  the  two  layers),  which  must  be  much  less  than 
1  in  order  for  linear  theory  to  be  applicable.  For  Loch  Ness 
in  September,  he  found  6  to  be  126a,  with  a  (the  half  amplitude 
of  the  seiche  motion  whose  wavelength  is  twice  the  lake  length) 
measured  in  meters.  He  thus  concluded  that  a  non-linear  ap¬ 
proach  must  be  used  to  explain  the  seiches  observed. 

A  qualitative  examination  of  the  temperature  data  shown 
in  Figures  9  and  10  indicates  that  the  seiche  period  internal 
waves  seen  in  Seneca  Lake  are  non-linear.  Amplitudes  of  over 
25  feet  are  seen  near  the  top  of  the  thermocline  when  the 
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epilimnion  is  only  about  50  feet  thick.  These  waves  cannot 
be  considered  small  amplitude,  so  simple  linear  theory  does 
not  apply.  To  show  this  quantitatively,  the  expansion  para¬ 
meter,  6,  was  calculated  for  the  typical  stratification  found 
in  Seneca  Lake  in  mid-July  and  mid-October.  The  July  value 
was  found  to  be  over  1000a  and  the  October  value  about  500a, 
with  a  measured  in  feet.  Since  the  amplitudes  seen  in  Seneca 
Lake  are  usually  well  over  ten  feet,  the  expansion  parameter 
will  be  much  greater  than  unity  and  non-linear  theory  must  be 
used. 

Thorpe  calls  these  non-linear  seiches  internal  undular 

bores  or  surges.  It  is  suggested  that  the  term  internal  surges 

be  used  to  describe  the  non-linear  seiches  as  the  term  bore  is 

usually  associated  with  tides.  Thorpe  gives  a  solution  for 

the  case  of  a  simple  harmonic  body  force  applied  to  the  upper 

p  h  2 

layer.  He  states  that  if  parameter  y=  —  - — ly  is  less  than  1, 

P2h  i 

the  surge  will  be  a  rise  in  the  level  of  the  interface 

and  if  y  is  greater  than  1  it  will  be  a  fall  in  the  level. 
Throughout  most  of  the  summer  and  fall,  y  will  be  greater  than 
1  in  Seneca  Lake  and  the  records  do  show  a  fall  in  the  level 
of  the  interface  (this  is  just  what  the  sharp  temperature  rises 
seen  in  Figure  9  mean;  a  moving  disturbance  with  a  sharp  front 
in  which  the  thermocline  is  lower  behind  the  front  will  cause 
a  stationary  thermistor  to  show  a  sharp  temperature  rise  as 
the  level  at  which  the  thermistor  is  changes  from  being  below 
the  thermocline  to  being  above  it  as  the  front  goes  by) . 

Thorpe  (1968A,  1968B  and  1971)  did  model  experiments  using 
a  rectangular  tank  filled  with  two  fluids  of  different  densities. 
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By  rocking  the  tank  (modelling  a  harmonic  body  force)  or 
tilting  it  and  then  bringing  it  to  rest  horizontally  (model¬ 
ling  the  relaxation  of  a  wind  stress)  he  was  able  to  produce 
surges . 

A  long,  narrow  plastic  rectangular  tank  with  inner  dimen¬ 
sions  of  6.5,  15.25  and  60.0  cm.  was  constructed  and  filled 
with  alcohol  and  turpentine.  The  tank  was  tilted,  causing 
the  interface  to  be  at  an  angle  with  the  top  and  bottom  of  the 
tank.  This  might  be  the  situation  in  a  long  narrow  lake  after 
a  wind  had  been  blowing  parallel  to  the  length  of  the  lake  for 
some  time.  The  tank  was  then  restored  to  the  horizontal  as  an 
analogy  to  the  lake  wind  dying  off  or  reversing  direction. 

The  resulting  surge  can  be  seen  in  Figure  17  in  which  it  is 
travelling  to  the  left.  In  the  figure  a  train  of  waves  can 
also  be  seen  to  be  travelling  with  the  surge  front  although 
these  might  not  be  analogous  to  wave  trains  seen  in  Seneca  Lake. 
In  the  tank  the  waves  have  sharp  crests  and  broad  troughs  while 
in  the  lake  the  reverse  is  seen.  Thorpe  (1971)  found  similar 
wave  trains  associated  with  surges  produced  in  his  tank  experi¬ 
ments.  He  was  unable  definitely  to  identify  similar  wave 
trains  in  his  records  of  surges  in  Loch  Ness. 

The  temperature  records  from  Seneca  Lake  clearly  show 
these  wave  trains;  indeed  their  amplitudes  in  many  cases  are 
greater  than  that  of  the  surge  itself  (see  Figures  18  through 


52)  . 
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shadow  of  the  true  interface  against  the 
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SHORT- PERIOD  INTERNAL  WAVES 

Introduction 

The  wave  train  seen  travelling  with  the  surge  in  the  tank 
model  in  Figure  17  and  similar  wave  trains  seen  by  Thorpe 
(197±)  in  his  model  experiments  are  very  predominant  in  Seneca 
Lake.  The  equations  governing  similar  wave  trains  associated 
with  non-linear  surface  bores  have  been  investigated  by  a  num¬ 
ber  of  authors.  Theoretical  profiles  of  these  waves  are  given 
by  Peregrine  (196*6) ,  Kadomstev  and  Karpman  (1971)  and  Witting 
(1972)  who  also  showed  photographs  of  model  experiments. 

The  large- amplitude  tascillations  with  periods  of  tens  of 
minutes  seen  in  Seneca  Lake  are  always  associated  with  the  in¬ 
ternal  surges  seen  in  the  lake.  The  occur  only  after  the  water 
temperatures  have  fallen  to  relatively  low  values  (i.e.,  the 
thermocline  has  risen)  and  after  they  pass,  water  temperatures 
are  higher  (the  thermocline  has  shifted  downward) .  This  can  be 
seen  in  Figure  10  in  which  the  temperature  data  obtained  during 
October  1968  and  digitized  at  one-hour  intervals  is  shown. 

Four  periods  of  large- amplitude  oscillations  were  seen  that 
month  and  the  initial  times  are  shown  by  arrows  in  the  figure. 
Figures  18  through  22  show  short-period  oscillations  recorded 
in  1969.  It  can  be  seen  from  these  figures  that  although  the 
short-period  oscillations  occur  only  in  association  with  a  long- 
period  surge,  not  all  of  the  surges  have  short-period  wave 
trains  associated  with  them  (e.g.,  the  surge  seen  at  about  470 
hours  in  Figure  10  and  the  one  seen  at  about  5462  hours  in 
Figure  18.  The  latter  has  one  sharp  spike  at  50  feet;  it  may 
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site  1  from  August  15,  1969,  1200  EST  to  August  19,  1969,  1200  EST. 
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Figure  19.  Similar  to  Figure  18.  The  data  are  from  September  30,  1969,  0000  EST 

to  October  6,  1969,  0000  EST. 
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1969,  0000 
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be  that  the  surge  was  not  developed  strongly  enough  to  generate 
a  complete  wave  train)  . 

Survey  of  the  short-period  oscillations  seen  in  Seneca  Lake 

Short-period  oscillations  were  seen  in  all  four  years  that 
thermistors  were  installed  on  the  lake.  The  earliest  recorded 
oscillations  that  could  be  identified  with  a  surge  occurred  in 
early  July,  1970.  The  oscillations  were  observed  throughout 
the  rest  of  each  summer  and  early  fall,  occurring  most  frequen¬ 
tly  in  October.  In  November,  as  the  temperature  differences 
between  the  epilimnion  and  the  hypolimnion  decreases,  the  sur¬ 
ges  and  the  associated  wave  trains  become  progressively  more 
indistinct. 

Figure  23  shows  the  surge  and  wave  train  of  July  15,  1971. 
The  data  from  the  two  25-foot  thermistors  show  large  tempera¬ 
ture  differences  due  to  instrumental  effects  (see  the  section 
on  temperature  probes).  Part  of  the  temperature  differences 
can  also  be  due  to  discrepancies  in  the  depth  of  the  thermistor 
pairs.  Because  it  shows  more  temperature  fluctuations  it  is 
probable  that  the  25-foot  thermistor  on  the  south  string  was 
actually  somewhat  deeper  than  its  counterpart  on  the  north 
string.  This  would  partly  explain  the  lower  temperature  recor¬ 
ded  by  that  probe.  Large- amplitude  oscillations  are  seen  only 
at  the  50-foot  level.  The  25-  and  75-foot  levels  show  some 
activity  while  the  lower  levels  show  almost  no  temperature  flue 
tuations .  This  does  not  necessarilv  mean  that  there  was  no 
internal  wave  activity  at  the  lower  levels;  the  water  was  al¬ 
most  isothermal  below  75  feet  so  it  would  take  very  large  ampli 
tude  waves  to  show  any  appreciable  temperature  changes.  (This 
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characteristic  limitation  of  thermistors  to  indicate  much  in¬ 
ternal  wave  activity  in  nearly  isothermal  water  is  graphically 
shown  in  Figure  24.  In  this  figure,  the  response  of  thermis¬ 
tors  at  five  depths  (25,  50,  75,  100  and  150  feet)  to  sinusoi¬ 
dal  internal  waves  with  amplitudes  linearly  increasing  from  0 
to  25  feet  peak-to-peak  are  shown  for  four  of  the  temperature 
structures  plotted  in  Figure  3.  The  sinusoidal  internal  waves 
used  to  construct  this  figure  had  the  same  amplitude  through¬ 
out  the  water  mass.  This  is  not  the  case  in  the  lake  where 
the  maximum  amplitude  will  be  in  the  thermocline  and  amplitudes 
will  decrease  both  above  and  below  it.  The  response  curves  for 
the  structure  of  September  25,  1969  show  flat-topped  tempera¬ 
ture  peaks  at  the  50-foot  level.  This  is  caused  by  the  (hypo¬ 
thetical,  in  this  case)  thermistor  entering  the  upper  mixed 
layer  and  is  commonly  seen  on  the  records.)  There  are  some 
small-amplitude  oscillations  seen  before  the  surge  front  passes 
the  thermistor  chains.  The  large-amplitude  oscillations  are 
initially  very  coherent  between  the  two  thermistor  strings  (see 
Table  1)  and  become  progressively  less  so  further  from  the 
surge  front,  indicating  increasing  turbulence.  This  was  the 
general-  pattern  seen  for  all  of  the  1971  surges  analyzed  for 
coherency  and  presumably  is  the  pattern  for  all  the  surges. 

The  amplitude  at  the  50-foot  level  is  about  25  feet,  peak-to- 
peak.  The  period  of  the  initial  large  oscillations  was  about 
9 . 4  minutes . 

Figure  25  shows  the  surge  front  of  July  28,  1971.  There 
is  more  activity  at  25  and  75  feet  than  during  the  previous 
surge.  The  beginning  of  the  wave  train  is  not  as  distinct  as 
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TABLE  1 :  COHERENCY  BETWEEN  NORTH  AND  SOUTH  THERMISTORS 
IN  1971 


Da  te 

Time  (EST) 

50 

Depth  in 
75 

feet 

100 

150 

Jul 

15 

1700-1750 

0.956 

1800-1900 

0.889 

1700-1900 

0.919 

1900-0000* 

0.778 

Jul 

28 

1200-1500 

0.805 

0.898 

1500-2055 

0.553 

0.908 

Jul 

30 

2300-2355 

0.937 

0.964 

Aug 

18 

1715-2035 

0.923 

0.924 

2030-0415* 

0.745 

0.630 

Sep 

22 

0400-0730 

0.991 

0.960 

0.958 

0730-1235 

0.919 

0.817 

0 . 874 

Sep 

24 

1210-1800 

0.783 

0 . 866 

1800-2100 

0.952 

0.950 

0.935 

2100-0140* 

0.944 

0.865 

0.779 

Oct 

3 

1500-1800 

0.976 

0.982 

1800-2340 

0.919 

0.799 

Oct 

6 

0800-1100 

0.985 

0.980 

0.975 

1100-1400 

0.840 

0.907 

0.897 

1400-1835 

0.743 

0.838 

Oct 

23 

1540-1650 

0.991 

0.986 

1650-1800 

0.978 

0.987 

0.975 

1500-1900 

0.989 

0.987 

0.952 

1900-0250* 

0.904 

0.921 

Oct 

27 

1300-1700 

0.958 

0.976 

0.974 

0.916 

1700-2020 

0 . 801 

0.914 

0.911 

Nov 

12 

1620-2020 

0.953 

0.989 

0.987 

2000-2350 

0.795 

0.948 

Nov 

22 

0210-0700 

0.456 

0700-0900 

0.928 

0900-1550 

0.847 

*of  the  following  day 
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in  most  cases.  The  coherency  at  the  50-foot  level  is  the 
lowest  for  any  of  the  1971  surges  analyzed. 

Figure  26  shows  the  surge  of  July  30,  1971.  The  drive 
cord  on  the  recording  instrument  broke  soon  after  the  front 
passed  so  only  a  portion  of  the  wave  train  was  recorded.  There 
was  considerable  activity  at  the  75-foot  level  in  addition  to 
that  at  the  50-foot  level.  The  25-  and  100-foot  levels  show 
some  oscillations  similar  to  that  seen  on  July  15,  1971  before 
the  front  arrives.  The  period  of  the  initial  large  oscil¬ 
lations  was  about  11.3  minutes. 

The  surge  of  August  18,  1971  is  shown  in  Figure  27.  This 
one  is  somewhat  different  from  those  seen  in  July.  There  is 
very  little  activity  before  the  front  arrives  and  it  contains 
three  distinct  oscillations  seen  from  25  to  100  feet.  Their 
period  was  about  13.0  minutes. 

The  surge  of  September  7,  1968  is  shown  in  Figure  30. 

There  is  some  small- amplitude  oscillation  seen  at  50  feet  be¬ 
fore  the  surge  arrives  but  the  front  is  sharp.  There  are  five 
temperature  peaks  seen  from  50  down  to  150  feet  with  the  ampli¬ 
tude  increasing  with  each  peak.  The  base  of  the  mixed  layer 
is  coming  down  below  50  feet  during  these  oscillations;  this 
is  the  explanation  of  the  flat- topped  peaks  seen.  The  maximum 
amplitude  of  these  oscillations  is  about  50  feet  and.  the  period 
is  about  15.9  minutes. 

Figure  31  shows  the  surge  of  September  22,  1971.  The  new¬ 
er  thermistor  strings  were  in  use  at  this  time  and  it  can  be 
seen  that  the  data  from  pairs  of  thermistors  at  the  same  depth 
appear  very  similar.  The  wave  train  of  this  surge  is  very 
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The  left  margin  at  5059  hours  is  July  30,  1971,  1900  EST. 
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The  left  margin  at  5508  hours  is  August  18,  1971.  1200  EST 
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1969 ,  1700  EST. 
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depths  are  25,  50,  75,  100  and  150  feet 
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Figure  30.  Similar  to  Figure  29  for  September  7,  1968. 
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Figure  31.  Similar  to  Figure  23.  The  left  margin  at  6336  hours  is  September  22,  1971 

0000  EST. 
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subdued.  There  is  a  peak  seen  at  50  feet  following  some  low- 
amplitude  activity.  The  maximum  that  brings  the  50-foot  level 
into  the  mixed  layer  is  seen  at  75  and  100  feet  but  there  is 
not  much  activity  thereafter  for  a  few  hours.  The  amplitude 
of  this  oscillation  is  about  30  feet  at  the  75-foot  level. 

The  surge  of  September  24,  1971,  shown  in  Figure  32,  is 
unusual,  there  being  much  activity  at  the  50-foot  level  before 
the  arrival  of  the  surge.  The  amplitudes  of  the  oscillations 
are  low,  water  motions  are  considerably  less  than  25  feet. 

The  surge  of  October  2,  1969  is  shown  in  Figure  33.  Both 
the  25-  and  50-foot  levels  are  seen  to  be  in  the  mixed  layer. 
The  75-foot  level  experienced  temperature  oscillations  of  al¬ 
most  10 °C  but  these  represent  water  motions  of  only  about  25- 
foot  amplitude.  Commencing  about  1830  EST  and  lasting  for 
over  two  hours,  oscillations  with  periods  less  than  5  minutes 
are  seen  at  the  75-foot  level.  (The  period  of  the  Brunt- 
Vaisala  frequency  at  this  level  is  also  just  under  5  minutes.) 
The  period  of  the  initial  oscillations  associated  with  the 
surge  front  is  about  16.4  minutes. 

Figure  36  shows  the  surge  of  October  6,  1971.  The  ampli¬ 
tude  is  about  30  feet.  The  50-foot  level  remains  out  of  the 
mixed  layer  until  a  number  of  hours  after  the  surge  passes. 

This  surge  contains  11  fairly  regular  initial  oscillations,  the 
most  seen.  The  period  is  about  17.6  minutes. 

The  surge  of  October  8,  1968  is  seen  in  Figure  37.  In 
many  respects  it  is  very  similar  to  that  of  October  4,  1969. 
Both  show  vertical  water  motions  of  over  50  feet  in  amplitude 
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Figure  36.  Similar  to  Figure  23.  The  left  margin  at  6674  hours  is  October  6,  1971, 
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Figure  37.  Similar  to  Figure  29  for  October  8,  1968. 
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at  the  100-foot  level  and  in  both  surges  the  maximum  amplitude 
is  seen  in  the  third  temperature  peak.  This  surge  also  shows 
regular  oscillations  at  the  75-foot  level  beginning  about  150 
minutes  after  the  front  passes  and  continuing  for  about  90 
minutes.  The  period  of  the  initial  oscillations  is  about  13.6 
minutes  and  that  of  the  later,  regular  ones  varies  from  about 
10  to  15  minutes. 

The  surge  of  October  13,  1969,  shown  in  Figure  40,  con¬ 
tains  regular  oscillations  beginning  about  4-1/2  hours  after 
the  front  came  through  and  lasting  more  than  1-1/2  hours,  simi¬ 
lar  to  that  of  October  8,  1968.  The  amplitude  of  the  initial 
oscillations  is  only  about  25  feet  and  the  period  about  22.5 
minutes . 

The  surge  of  October  18,  1969,  seen  in  Figure  42,  is  some¬ 
what  unusual  in  appearance.  There  are  short-period  oscilla¬ 
tions  at  75  feet  similar  to  that  seen  on  September  24,  1971. 

The  150-foot  thermistor  indicates  a  wave  amplitude  of  over  50 
feet  at  that  depth,  yet  the  75-foot  level  is  below  the  mixed 
layer  throughout  most  of  this  record.  The  period  of  the  ini¬ 
tial  oscillations  is  about  19.7  minutes. 

The  surge  of  October  22,  1969,  seen  in  Figure  43,  begins 
differently  from  the  previous  ones  seen  in  October.  The 
initial  temperature  rise  is  not  as  rapid  and  the  first  decline 
not  as  smooth  nor  as  rapid  as  in  previous  surges.  The  75- foot 
thermistor  is  probably  above  the  base  of  the  mixed  layer  for 
most  of  the  record  although  from  the  figure  it  appears  to  be 
at  a  lower  temperature  than  the  25-  and  50- foot  ones.  This  is 
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Figure  39.  Similar  to  Figure  29  for  October  12,  1968. 
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Figure  40.  Similar  to  Figure  28.  The  left  margin  at  6845  hours  is  October  13,  1969, 
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probably  due  to  the  calibration  difficulties  experienced  with 
the  older  thermistor  chain. 

Figure  47  shows  the  surge  of  October  30,  1968.  Disregar¬ 
ding  the  fact  that  temperature  peaks  at  the  100- foot  level  are 
higher  than  even  the  25-foot  temperatures  (this  is  due  to  the 
calibration  difficulties  encountered  with  the  old  thermistor 
string) ,  it  is  apparent  that  the  flat  temperature  peaks  at 
this  level  indicate  that  the  mixed  layer  is  oscillating  to 
below  100  feet.  This  surge  is  unusual  for  the  oscillations 
seen  at  75  feet  before  the  front  arrives  and  because  of  the 
single  temperature  peak  seen  as  the  surge  arrives. 

The  surge  of  November  4,  1969  is  shown  in  Figure  48.  The 
mixed  layer  is  below  75  feet  before  the  arrival  of  the  surge 
and  oscillates  above  and  below  100  feet  after  the  front  comes 
through.  Oscillations  are  close  to  50  feet  in  amplitude  at 
the  150- foot  level  and  the  period  is  the  longest  observed, 
about  40.3  minutes. 

The  oscillations  associated  with  the  surge  front  of  Novem 
ber  16,  1968,  shown  in  Figure  52,  build  up  slowly  in  contrast 
to  those  seen  earlier  in  the  year.  Fluctuations  are  seen  at 
50  feet  which  is  very  unusual  for  so  late  in  the  year. 

Figure  53  shows  the  temperature  data  of  November  22,  1971 
It  is  questionable  as  to  whether  this  represents  a  surge. 

It  is  evident  from  the  figures  that,  although  exceptions 
can  be  found  for  almost  any  generalization  made  about  the  sur¬ 
ges,  there  are  a  number  of  characteristics  that  are  common  to 
all  or  most  of  the  surges  seen  in  Seneca  Lake.  In  all  cases, 
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string.  The  left  margin  at  7090  hours  is  October  23,  1971,  1000  EST 
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Figure  45.  Similar  to  Figure  29  for  October  25,  1968. 
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Figure  47.  Similar  to  Figure  29  for  October  30,  1968. 
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Figure  48.  Similar  to  Figure  28.  The  left  margin  at  7378  hours  is  November  4,  1969, 
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igure  49.  Similar  to  Figure  29  for  November  7,  1968. 
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Figure  51.  Similar  to  Figure  44.  The  left  margin  at  7570  hours  is  November  12,  1971 
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Figure  52.  Similar  to  Figure  29  for  November  16,  1968. 
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temperatures  have  been  slowly  dropping  at  all  levels  below 
the  mixed  layer  (the  temperature  structure  curve  is  shifting 
upward)  for  tens  of  hours  before  the  surge  is  seen,  indi¬ 
cating  that  the  warm  water  of  the  epilimnion  is  moving  away 
and  is  being  replaced  by  colder  water  from  below.  In  all 
cases,  the  first  strong  oscillation  is  a  rise  in  temperature 
and  when  the  wave  train  had  passed,  all  levels  showed  either 
the  same  or  a  higher  temperature  than  before  the  surge.  This 
shows  that  the  surge  is  always  a  fall  in  the  level  of  iso¬ 
therms,  indicating  warm  water  moving  in  and  colder  water  being 
displaced  downward.  In  most  cases,  the  first  oscillations 
start  rather  suddenly  and  the  first  few  are  smooth  and  regular. 
In  all  cases,  the  frequency  of  the  initial  oscillations  are 
less  than  the  Brunt- Vaisala  frequency.  The  first  few  waves 
generally  have  the  largest  amplitudes  and  are  very  coherent 
over  at  least  a  few  hundred  feet;  the  following  oscillations, 
which  usually  continue  for  hours  and  generally  are  irregular, 
are  less  coherent  and  contain  more  turbulence.  The  wave  trains 
seen  in  1971  were  generally  of  smaller  amplitude  than  those 
of  1968  and  1969  and  their  initial  temperature  rise  was  not  as 
sudden.  This  may  be  due  to  the  different  location  used  in  1971. 
Table  2  contains  a  summary  of  the  data  for  all  of  the  short- 
period  wave  trains  associated  with  surges  seen  in  the  years 
1968,  1969  and  1971. 

In  Figure  54  the  period  of  the  initial  oscillations  of 
the  wave  train  are  plotted  against  the  date  of  the  surge.  It 
can  be  seen  that  there  is  a  definite  increase  in  period  from 
July  through  November.  A  straight  line  can  be  drawn  through 
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TABLE  2:  SUMMARY  OF  DATA  FOR  SHORT-PERIOD  INTERNAL  WAVES 


Date 

Initial 

period 

(min) 

Mixed  layer  depth 
before  after 

(ft)  (ft) 

Max 
ampl . 
(ft) 

Depth  of 
max  ampl 
(ft) 

J  ul 

15/71 

9 . 4  ±0 . 5 

<25 

<25 

25 

50 

Jul 

28/71 

<25 

<25 

10 

50 

Jul 

30/71 

11. 3±0 . 5 

<25 

25 

50 

Aug 

18/71 

13.0±0.5 

<25 

<25 

20 

50 

Aug 

18/69 

14.110.5 

25 

35 

15 

50 

Aug 

24/68 

17.810.5 

30 

40 

20 

50 

Sep 

7/6  8 

15.9+0.5 

25 

35 

50 

75 

Sep 

22/71 

30 

45 

30 

75 

Sep 

24/71 

30 

40 

10 

75 

Oct 

2/69 

16.410.5 

55 

70 

25 

75 

Oct 

3/71 

30 

40 

20 

75 

Oct 

4/69 

11.710.5 

30 

50 

60 

100 

Oct 

6/71 

17.610.5 

30 

50 

30 

75 

Oct 

8/68 

13.610.5 

30 

65 

65 

100 

Oct 

9/69 

18.810.5 

45 

55 

25 

100 

Oct 

12/68 

16.910.5 

60 

70 

65 

100 

Oct 

13/69 

22.5+0.5 

40 

50 

25 

100 

Oct 

15/69 

19.210.5 

35 

50 

50 

100 

Oct 

18/69 

19.710.5 

50 

60 

60 

150 

Oct 

22/69 

50 

70 

40 

100 

Oct 

23/71 

21.110.5 

45 

60 

20 

100 

Oct 

25/68 

19.7+0.5 

50 

75 

50 

100 

Oct 

27/71 

35 

55 

20 

100 

Oct 

30/68 

60 

85 

25 

100 

Nov 

4/69 

40.310.5 

75 

100 

40 

150 

Nov 

Nov 

7/68 

10/69 

31.111.5 

26.311.0 

80 

90 

30 

100 

Nov 

Nov 

12/71 

16/68 

30.012.0 

60 

80 

35 

150 

Nov 

22/71 

110 

120 

30 

150 
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Figure  54.  Period  of  the  initial  oscillations  of  the  wave  trains 

associated  with  internal  surges  seen  in  the  years 
1968 ,  1969  and  1971.  The  vertical  broken  lines  delin¬ 
eate  separate  months  beginning  with  July. 
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the  points  from  July  through  October  and  it  would  show  an 
increase  in  period  of  1  minute  about  every  12  days.  However, 
there  would  be  a  number  of  points,  especially  in  October,  that 
would  not  fit  the  line  well.  The  periods  measured  in  Novem¬ 
ber  were  much  higher  than  in  previous  months  and  do  not  fit  a 
linear  pattern.  An  exponential  curve  of  initial  period  versus 
time  of  year  would  fit  the  total  data  better  but  there  are  not 
enough  points  to  attempt  a  determination  of  the  curve  para¬ 
meters.  The  initial  period  is  probably  a  function  of  the  ther¬ 
mal  structure  of  the  lake  which  also  changes  throughout  the 
summer  and  fall  (see  Figure  3)  although  no  simple  relationship 
between  the  various  parameters  (densities,  layer  thicknesses) 
and  the  period  has  been  found. 

Another  way  to  present  the  temperature  data  is  to  plot 
the  depths  of  several  isotherms  as  a  function  of  time  rather 
than  plotting  the  original  temperatures  at  specific  depths  ver¬ 
sus  time.  This  type  of  presentation  is  shown  for  four  of  the 
surges  seen  in  1968  in  Figures  55  through  58.  A  linear  inter¬ 
polation  was  used  to  determine  temperatures  between  thermistors. 
This  distorts  the  isotherms  drawn  in  the  figures  because  the 
real  temperature  structure  is  more  like  an  exponential  function 
of  depth  than  a  linear  one.  Also,  because  of  calibration  dif¬ 
ficulties,  comparisons  of  temperatures  between  two  thermistors 
at  different  depths  cannot  be  done  very  accurately.  Although 
only  illustrative,  Figures  55  through  58  show  more  clearly  that 
the  surges  are,  in  all  cases,  a  lowering  of  the  isotherm  levels. 
The  figures  also  show  that  the  short-period  oscillations  have 
sharp  troughs  and  flat  crests,  in  contrast  to  the  waves  seen 
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Figure  56.  Isotherm  depths  for  surge  of  September  7,  1968  (Figure  30). 
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Figure  58.  Isotherm  depths  for  surge  of  November  7,  1968  (Figure  49). 
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in  the  tank  model  (Figure  17) . 

Power  spectra  of  the  short-period  wave  trains  were  ob¬ 
tained  for  three  of  the  surges  seen  in  1969.  Twenty  hours  of 
data  at  one-minute  intervals,  beginning  shortly  before  the 
surge  was  seen,  were  used  to  obtain  the  spectra.  Figure  59 
shows  the  power  spectrum  of  temperature  at  50  feet  for  the 
surge  of  August  18,  1969  seen  in  Figure  28.  There  is  a  spec¬ 
tral  peak  at  15.0  minutes  which  is  close  to  the  period  mea¬ 
sured  from  the  records,  14.1  minutes.  The  power  spectrum  of 
October  13,  1969  at  75  feet  is  shown  in  Figure  60.  There  are 
two  peaks  seen  on  this  spectrum;  the  first  at  26.67  minutes 
is  close  to  that  measured  for  the  initial  oscillations,  while 
the  second  peak  at  10.91  minutes  corresponds  to  the  period  of 
the  oscillations  seen  later  during  this  surge  (see  Figure  40) . 
The  power  spectrum  at  150  feet  of  the  surge  of  November  4, 

1969  (Figure  48)  is  shown  in  Figure  61.  There  is  a  peak  at 
40.0  minutes,  which  corresponds  well  with  the  measured  period 
of  40.3  minutes.  In  all  three  spectra  the  peaks  discussed  are 
not  above  the  background  by  much  more  than  the  90%  confidence 
band. 

The  two  thermistor  chains  used  in  1971  at  site  2  were 
used  to  determine  the  direction  of  travel  of  the  surges.  In 
all  12  surges  seen  that  year,  the  initial  oscillations  were 
always  seen  at  the  south  end  of  the  SMP  before  they  were  seen 
at  the  north  end  (see  Figures  23,  25,  26,  27,  31,  32,  34,  36, 
44,  46,  51  and  53).  The  spacing  of  the  two  thermistor  strings 
on  the  SMP,  however,  was  not  great  enough  to  determine  phase 
velocities  accurately.  Phase  velocities  were  determined  for 
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Figure  59.  Power  spectrum  of  20  hours  of  data  at  one  minute  inter¬ 
vals  beginning  August  18,  1969,  1810  EST.  Thermistor 
depth  is  50  feet. 
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Figure  60.  Power  spectrum  of  20  hours  of  data  at  one  minute  inter¬ 
vals  beginning  October  13,  1969,  0640  EST.  Thermistor 
depth  is  75  feet. 
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Figure  61.  Power  spectrum  of  20  hours  of  data  at  one  minute  inter¬ 
vals  beginning  November  4,  1969,  1020  EST.  Thermistor 

depth  is  150  feet. 


103 


the  five  surges  that  were  seen  at  both  sites  1  and  2,  and  are 
listed  in  Table  3.  The  data  for  the  surge  of  October  21,  1971 
are  puzzling  at  first  glance  because  the  arrival  times  of  the 
first  temperature  peak  give  phase  velocities  considerably 
higher  than  those  determined  from  the  second  and  third  peaks. 
Inspection  of  Figure  45  shows,  however,  that  the  first  few  oscil¬ 
lations  of  this  surge  are  somewhat  unusual.  There  are  at  least 
three  possibilities:  (1)  A  temperature  peak  is  missing  between 
the  first  two  seen  at  the  SMP.  If  this  is  the  case,  then  the 
second  peak  seen  at  the  SMP  should  be  compared  with  the  third 
at  the  TCP,  yielding  a  travel  time  of  2  hours,  11  minutes  and 
a  phase  velocity  of  49.24  cm. /sec.  (2)  The  first  peak  seen  at 
the  SMP  has  emerged  from  the  surge  front  during  the  time  of 
travel  between  the  TCP  and  SMP.  If  this  is  the  case,  then  the 
second  SMP  peak  should  be  compared  with  the  first  one  seen  at 
the  TCP  and  the  third  SMP  peak  with  the  second  TCP  peak.  This 
yields  travel  times  ranging  from  2  hours ,  51  minutes  to  2  hours , 
57  minutes  and  phase  velocities  ranging  from  36.44  to  37.72 
cm. /sec.  (3)  A  combination  of  (1)  and  (2)  above,  i.e.,  a  tem¬ 
perature  peak  is  missing  on  the  SMP  and  the  first  peak  seen 
there  has  emerged  after  the  surge  passed  the  TCP.  If  this  is 
the  case,  then  the  second  SMP  peak  should  be  compared  with  the 
second  TCP  peak  and  the  third  with  the  third.  This  is  listed 
in  Table  3,  the  resulting  phase  velocities  are  between  42.16 
and  43.25  cm. /sec. 

A  comparison  of  the  first  three  peaks  seen  on  the  records 
obtained  at  site  1  with  those  of  site  2  suggests  that  the 
third  hypothesis  given  above  is  correct,  however,  the  compari- 


TABLE  3:  TRAVEL  TIMES  OF  SURGES  BETWEEN  SITES 
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son  was  not  conclusive  because  the  chart  paper  at  the  TCP  ran 
out  at  the  time  of  arrival  of  the  fourth  peak. 

The  observed  speeds  determined  for  the  three  other  surges 
in  which  more  than  the  first  temperature  maximum  was  used  are 
internally  consistent  to  within  1  cm. /sec.  Timing  discrepan¬ 
cies  between  the  two  barges  are  known  to  be  less  than  five 
minutes  at  the  worst,  so  it  is  felt  that  the  measured  phase 
velocities  listed  in  Table  3  are  certainly  accurate  to  better 
than  ±2  cm/sec.  In  the  last  column  of  Table  3  are  listed  the 
theoretical  travel  times  of  the  surges,  calculated  using  linear 
theory  and  a  multilayered  lake  model  (see  Appendix  II) .  It 
is  seen  that  the  calculated  and  observed  speeds  are  very  close 
for  the  three  cases  (July  15,  October  23  and  November  12)  in 
which  the  measured  speeds  are  considered  very  reliable,  but 
there  does  not  seem  to  be  a  simple  relationship  between  the 
two.  In  particular,  the  Froude  number,  which  is  the  ratio  of 
the  surge  speed  to  that  of  a  linear  wave  and  which  theoreti¬ 
cally  is  always  greater  than  one,  would  be  less  than  unity  for 
two  of  those  surges.  A  possible  explanation  of  this  discre¬ 
pancy  will  be  attempted  in  the  following  section. 

Theoretical  and  experimental  studies  of  non-linear  waves 

It  has  already  been  shown  that  the  long-period  internal 
waves  seen  in  Figures  9  through  22  are  non-linear  surges  and 
cannot  be  understood  in  terms  of  simple  linear  models.  What 
of  the  short-period  wave  trains  seen  in  Figures  23  through  53? 
Many  of  the  initial  oscillations  seen,  while  not  sinosoidal 
in  shape,  are  symmetrical  with  respect  to  time.  Many  also 
have  a  period  that  can  be  measured  (Table  2) .  A  first  attempt 
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was  made  to  understand  the  wave  trains  as  linear  waves  gener¬ 
ated  by  the  non-linear  surge  and  coupled  to  it  by  the  phase 
velocity.  A  matrix  technique  similar  to  that  employed  by 
Haskell  (1951)  for  Rayleigh  and  Love  waves  was  used  to  obtain 
dispersion  curves  for  internal  waves  in  many-layered  fluids. 

In  the  course  of  the  investigation,  however,  it  became  evident 
that  the  wave  trains  associated  with  the  surges  are  non-linear 
and  the  matrix  method  developed  did  not  apply.  Because  this 
method  can  be  useful  in  other  problems  of  internal  waves,  a 
summary  of  the  technique  is  given  in  Appendix  II. 

As  pointed  out  earlier  in  this  report,  a  number  of  authors 
have  seen  short-period  wave  trains  travelling  with  internal 
surges  generated  in  tank  models.  Figure  17  is  a  photograph  of 
one  case.  Several  authors  have  theoretically  investigated  the 
evolution  of  an  initial  perturbation  on  the  surface  of  a  fluid 
using  the  Korteweg-de  Vries  equation  (Korteweg  and  de  Vries , 
1895)  and  have  shown  that  wave  trains  similar  to  those  seen  in 
laboratory  experiments  and  in  Seneca  Lake  develop  from  a  surge 
in  certain  circumstances.  Peregrine  (1966)  used  a  numerical 
technique  to  investigate  the  equations  and  shows  figures  illus¬ 
trating  the  results.  In  particular,  he  shows  the  growth  and 
development  of  undular  bores  under  several  initial  conditions 
(see  his  Figures  3  through  6  which  look  very  much  like  the 
surges  seen  in  Seneca  Lake) . 

Witting  (1972)  classified  five  types  of  surface  bores. 

He  shows  photographs  and  illustrative  figures  of  several  types. 
The  different  types  of  bores  occupy  separate  areas  on  a  graph 
of  bore  strength  (which  is  measured  by  the  Froude  number,  F^, 
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which  is  the  ratio  of  the  bore  speed  to  that  of  a  linear  wave) 
versus  the  depth  of  w ater.  A  Froude  number  can  be  calculated 
for  the  three  surges  listed  in  Table  3  whose  measured  speeds 
are  considered  reliable,  but  two  of  them  would  be  less  than 
one.  A  possible  explanation  of  the  situation  is  suggested  by 
the  positions  of  the  two  observation  sites  in  Figure  1.  It  is 
seen  from  the  map  that  the  sites  are  near  a  bend  in  the  lake. 

A  straight-line  distance  between  the  sites  was  used  to  obtain 
the  phase  velocities  listed  in  Table  3.  It  is  probable,  how¬ 
ever,  that  a  surge  does  not  follow  a  straight  line  ray  path 
between  the  two  barges.  This  is  shown  schematically  in  Figure 
62.  In  the  upper  portion  of  the  figure  is  shown  the  case  in 
which  the  two  barges  are  not  on  the  same  ray  path,  which  is 
assumed  to  be  straight.  In  this  case,  the  true  distance  tra¬ 
velled  by  a  wave  front  would  be  less  than  that  measured  by  a 
factor  of  the  cosine  of  the  angle  between  the  ray  path  and  the 
direct  path  between  the  barges.  This  would  reduce  the  measured 
speeds  and  worsen  the  situation.  The  lower  portion  of  the 
figure  shows  the  case  in  which  the  two  barges  are  on  the  same 
ray  path  but  it  is  curved  (in  the  figure  it  is  a  75  degree 
circular  arc,  which  has  enough  curvature  to  result  in  Froude 
numbers  greater  than  one  for  the  three  surges  being  considered) . 
In  this  case,  the  distance  between  the  barges  along  the  ray 
path  would  be  greater  than  the  straight-line  distance  and  the 
measured  speeds  would  be  increased.  The  actual  ray  path  of 
a  surge  between  the  two  sites  might  be  a  combination  of  both 
cases.  If  the  second  predominates,  the  problem  of  Froude  num¬ 
bers  apparently  less  than  unity  would  be  resolved.  It  is 
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Figure  62. 


Two  possible  ray  paths  between  the  barge  sites. 
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apparent  though,  that  with  the  data  obtained  very  little  can 
be  said  about  Froude  numbers  other  than  they  are  close  to  one. 

Kadomtsev  and  Karpman  (1971)  discuss  the  evolution  of  an 
initial  perturbation  into  a  set  of  solitary  waves  or  solitons. 
Their  Figure  4,  which  was  obtained  by  numerically  solving  the 
Korteweg-de  Vries  equation,  looks  much  like  the  initial  large- 
amplitude  oscillations  seen  in  the  lake.  When  dissipation  is 
included,  the  Korteweg-de  Vries-Burgers  equation  results.  Sol¬ 
utions  of  this  equation,  in  addition  to  producing  an  initial 
train  of  large -amplitude  waves,  also  result  in  a  change  in 
level  behind  the  initial  oscillations  similar  to  the  depression 
of  isotherms  following  a  surge  seen  in  the  lake.  Although  the 
evidence  obtained  thus  far  is  only  qualitative,  it  appears 
likely  that  the  surges  seen  in  Seneca  Lake  can  be  understood 
in  terms  of  nonlinear  waves  governed  by  the  Korteweg-de  Vries 
equation . 

In  summary,  the  temperature  data  obtained  in  Seneca  Lake, 
when  filtered  and  sampled  at  periods  of  a  few  hours  or  more, 
look  similar  to  that  obtained  from  other  lakes  by  many  workers. 
The  oscillations  seen  were  explained  as  temperature  seiches 
of  the  thermo cline  and  this  hypothesis  seemed  to  explain  most 
of  the  important  features  seen  in  the  data.  By  1953,  with 
the  publication  of  Mortimer's  paper,  the  subject  of  long-period 
internal  waves  in  long  narrow  lakes  seemed  to  be  well  under¬ 
stood.  However,  in  1971,  Thorpe  showed  that  the  internal  seiche 
in  Loch  Ness  appeared  to  be  a  nonlinear  surge  rather  than  a 
simple  standing  wave  controlled  by  the  dimensions  of  the  lake. 
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The  data  obtained  from  Seneca  Lake  show  conclusively  that  in 
this  lake,  at  least,  there  are  large- amplitude  undular  surges 
that  appear  to  be  the  internal  analogues  of  the  surface  bores 
investigated  by  several  authors.  It  is  therefore  apparent 
that  a  whole  new  avenue  of  investigation  is  opening  up. 

Wind  measurements 

In  the  classic  description  of  internal  seiches  in  lakes, 
the  wind  plays  a  role  of  initially  inclining  the  isotherms 
which  then  respond  by  rocking  back  and  forth  at  the  natural 
period  of  the  lake.  Briefly  stated,  the  sequence  of  events 
is,  according  to  this  theory:  (1)  a  wind  blows  roughly  para¬ 
llel  to  the  long  axis  of  the  lake  for  a  sufficient  time  to 
(2)  blow  a  sizeable  amount  of  the  warm  surface  water  to  one 
end  of  the  lake  with  the  colder,  deeper  water  replacing  it  at 
the  other  end.  The  thermocline  is  therefore  tilted,  the  tilt 
maintained  by  the  opposition  of  the  pressure  gradient  which 
tends  to  restore  the  thermocline  to  the  horizontal  and  the 
wind  stress  which  tends  to  push  even  more  warm  water  down  the 
lake;  (3)  after  a  period  of  time,  the  wind  dies  down.  The 
tilted  thermocline  is  no  longer  a  stable  condition  and  (in 
most  cases)  it  oscillates  up  and  down  at  the  ends  of  the  lake 
with  a  node  near  the  center.  These  simple  harmonic  oscil¬ 
lations  die  out  owing  to  frictional  effects.  Although  the 
internal  waves  seen  in  Seneca  Lake  have  the  character  of  in¬ 
ternal  surges  travelling  up  the  lake  rather  than  a  linear  oscil¬ 
lations  of  the  thermocline,  wind  stress  would  still  appear  to 
be  the  most  likely  cause  of  the  initial  imbalance  that  starts 
the  surge  going. 
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In  order  to  test  the  theory,  winds  were  monitored  at 
site  1  (Figure  1)  for  15  days  in  August,  1970,  and  for  two 
periods  in  1971,  totalling  63  days.  Wind  speed  and  direction 
were  recorded  on  strip  charts  and  digitized  at  one-hour 
intervals.  No  surges  were  seen  during  the  time  wind  data 
was  collected  in  1970  but  the  wind  data  itself  is  interesting. 
Regular  oscillations,  with  a  period  of  24  hours  are  seen  es¬ 
pecially  in  the  north-south  component  (see  Figure  63) .  These 
are  very  prominent  between  the  11th  and  16th  (days  222  to 
228  in  the  figure)  of  August.  A  close  examination  of  the 
data  reveals  maximum  north  velocities  in  the  afternoons  (usu¬ 
ally  before  6  p.m. ),  maximum  south  velocities  in  the  mornings 
(usually  around  6  a.m. )  with  zero  north-south  velocities  be¬ 
fore  noon  and  midnight.  The  geography  of  the  region  offers  a 
reasonable  explanation  of  the  wind  pattern.  Seneca  Lake  is 
in  a  valley  oriented  north-south  and  is  between  the  mountains 
of  the  Northern  Appalachian  Plateau  to  the  south  and  the  plain 
to  the  south  of  Lake  Ontario  to  the  north.  In  the  afternoon, 
the  north  winds  are  coming  from  the  direction  of  the  plains 
and  are  blowing  up  the  valley.  In  the  early  morning,  the 
south  winds  are  coming  from  the  direction  of  the  mountains  and 
are  blowing  down  the  valley.  This  is  just  what  is  expected  in 
the  classic  mountain  and  valley  wind  (see,  for  example,  Defant, 
1951)  and  it  is  therefore  concluded  that  this  is  the  normal 
wind  pattern  of  the  lake  and  is  seen  whenever  it  is  not  ob¬ 
scured  by  "meteorological"  winds. 

Figure  64  shows  the  wind  data  for  25  days  beginning  June 
The  mountain- valley  wind  is  also  seen  at  times  in 
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TIME  IN  DAYS  FROM  JAN-  1,  1970  0000  E5T 


Figure  63.  Wind  speeds  recorded  at  barge  site  1  from  August  7 ,  1970 

through  August  21,  1970.  The  north-south  component  is 
drawn  in  solid  line,  the  east-west  component  in  dashed 


line.  Positive  values  are  north  and  east. 
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this  figure  (e.g. ,  from  day  189  through  191).  Only  one  surge 
was  seen  during  the  period  shown  in  Figure  64;  it  passed  site 
2  on  July  15  (day  195)  at  1712  EST.  Figure  65  shows  the  wind 
data  for  38  days  beginning  October  21.  Times  when  surges 
were  seen  are  marked  with  arrows;  solid  ones  for  those  seen 
at  both  sites  and  dashed  arrows  for  surges  seen  only  at  site 
1  (of  the  four  surges  seen  only  at  site  1,  three  were  probably 
seen  at  site  2  but  could  not  be  definitely  correlated  because 
difficulties  with  the  chart  drive  of  the  recorder  at  site  2 
made  it  impossible  to  determine  arrival  times  at  that  site) . 

It  is  apparent  that  a  simple  relationship  between  the 
wind  and  the  internal  surges  does  not  exist.  Some  surges  occur 
after  strong  wind  activity,  others  after  weak  activity,  while 
some  strong  winds  do  not  produce  surges  at  all. 

Another  possible  mechanism  for  piling  up  warm  water  at 
one  end  of  the  lake  is  the  atmospheric  pressure  gradient.  A 
pressure  difference  of  one  millibar  between  the  ends  of  the 
lake  would  produce  an  equilibrium  difference  in  surface  level 
of  about  one  centimeter.  Because  the  density  difference  be¬ 
tween  the  epilimnion  and  the  hypolimnion  is  on  the  order  of 
10  J  times  the  density  of  water,  a  difference  of  about  10 
meters  in  the  depth  of  the  thermocline  at  the  two  ends  of  the 
lake  could  be  produced.  While  the  usual  atmospheric  pressure 
differences  along  the  lake  is  less  than  one  millibar,  dif¬ 
ferences  of  that  magnitude  are  not  uncommon.  However,  as  with 
wind  speeds,  no  relationship  was  found  between  the  pressure 
gradient  and  surge  generation. 

It  is  interesting  that  all  the  surges  for  which  a  direc- 
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tion  of  travel  has  been  determined  were  found  to  be  moving 
north.  The  thermocline  is  in  continuous  motion  (see  Figures 
9  through  11  and  18  through  22) .  It  is  possible  that  the 
undular  surges  only  develop  at  the  southern  end  of  the  lake 
when  the  wind  re-enforces  a  weak  surge  already  travelling  the 
lake.  The  topography  of  the  lake  bottom  is  very  different  at 
the  two  ends  (Figure  66) .  The  southern  part  is  deeper  than 
the  northern  and  the  slope  is  much  steeper.  Less  than  0.6 
miles  from  the  south  end  of  the  lake,  the  water  is  over  200 
feet  deep  and  depths  over  400  feet  are  encountered  about  1.75 
miles  from  the  southern  end.  The  slope  is  much  gentler  at  the 
northern  end;  0.6  miles  from  shore  the  w ater  is  less  than  25 
feet  deep.  Water  depths  in  excess  of  200  feet  are  not  en¬ 
countered  until  about  3.7  miles  from  the  northern  end  and  400 
foot  depths  about  7.8  miles.  Thus,  a  weak  surge  travelling 
north  encounters  gradually  shoalling  water,  while  one  travel¬ 
ling  to  the  south  practically  hits  a  wall.  It  might  take  the 
proper  combination  of  the  strength  of  a  weak  surge,  its  reflec¬ 
tion  at  the  southern  end  of  the  lake,  re-enforced  by  proper 
winds  to  produce  the  large -amplitude  undular  surges  seen  near 


the  center  of  the  lake. 
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Figure  66.  Longitudinal  section  of  Seneca  Lake.  North  is  to  the  left.  Sites  1  and 

2  as  per  Figure  1.  Depths  are  in  feet,  distance  in  miles. 
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CONCLUSION 

Temperature  measurements  obtained  during  the  course  of 
four  years  have  shed  considerable  light  on  the  behavior  of 
Seneca  Lake.  Bathythermograph  records  have  been  used  to  re¬ 
construct  the  yearly  temperature  cycle.  These  records  indi¬ 
cate  that  the  lake  is  probably  freely  circulating  throughout 
much  of  the  winter  although  the  temperature  is  below  that  of 
maximum  density  water.  These  records  have  also  been  used  to 
calculate  heat  budgets.  The  heat  budget  curves  (Figures  4 
through  8)  indicate  that  as  the  upper  layers  of  the  lake  lose 
heat  in  autumn,  lower  layers  are  warming.  This  transfer  of 
heat  downward  may  by  accomplished  to  a  large  extent  by  the 
large- amplitude  internal  surges  seen  in  late  summer  and  early 
fall. 

These  large-amplitude  internal  waves  are  the  most  striking 
feature  of  the  temperature  data  seen.  Although  similar  in 
amplitude  to  the  internal  "temperature"  seiches  seen  in  long, 
narrow  lakes  since  the  beginning  of  this  century,  the  internal 
waves  seen  in  Seneca  Lake  are  non-linear  and  have  large-ampli¬ 
tude,  short-period  wave  trains  associated  with  them.  They 
closely  resemble  undular  surges  that  have  recently  been  inves¬ 
tigated  both  theoretically  and  with  tank  models.  It  is  con¬ 
cluded  that  the  internal  waves  seen  in  the  lake  are  internal 
analogues  of  these  surface  non-linear  waves  and  are  called 
internal  undular  surges. 

The  consequences  of  these  internal  surges  have  not  been 
investigated  but  they  could  be  of  major  importance.  There  is 
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probably  much  turbulent  motion  in  the  wave  train  associated 
with  a  surge.  This  turbulence  may  play  an  important  role  in 
the  downward  transport  of  heat  seen  in  late  summer  and  early 
autumn.  It  may  also  bring  nutrients  from  the  lower  layers 
into  the  upper  mixed  layer,  renew  oxygen,  etc.  The  consequen¬ 
ces  to  the  life  forms  in  the  lake  of  the  large  and  rapid 
oscillations  of  isothermal  surfaces  could  be  profound. 

The  wind  data  collected  in  1971  do  not  reveal  a  simple 
mechanism  for  the  generation  of  the  internal  surges.  Because 
all  of  the  surges  whose  direction  could  be  determined  were 
travelling  north,  bottom  topography  (which  is  very  different 
at  the  two  ends  of  the  lake)  must  play  an  important  role. 

The  non-linear  internal  undular  surges  seen  in  Seneca 
Lake  reopen  an  area  of  investigation  that  for  many  years  lim- 
nologists  felt  was  solved  except  for  minor  details.  Are  the 
surges  seen  in  this  lake  an  unusual  occurrence  or  are  they 
the  general  rule  for  narrow  lakes?  Thorpe  (1971)  has  iden¬ 
tified  non-linear  surges  in  Loch  Ness  but  has  not  seen  wave 
trains,  either  because  his  instruments  were  not  sensitive 
enough  or  because  they  do  not  exist  in  that  lake.  Other  lakes 
must  be  investigated. 

Work  should  also  continue  in  Seneca  Lake.  Because  of 
the  two  barge  sites  near  the  center  of  the  lake,  data  collec¬ 
tion  there  is  relatively  easy.  This  should  be  augmented  by 
data  collected  at  other  points  in  the  lake  so  that  surges 
could  be  followed  as  they  travel  through  it.  Investigators 
have  been  using  small  tanks  to  model  nature.  In  Seneca  Lake, 
nature  has  given  us  a  large  tank  in  a  convenient  location  to 
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play  with.  Advantage  should  be  taken  of  this. 


Appendix  I. 


Tables  of  Heat  Budget  Data 


TABLE  1.1:  HEAT  BUDGET  DATA  -  1968 
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Appendix  II.  Linear  Internal  Waves  in  Multilayered 
Fluids  -  A  Matrix  Technique  of  Obtain  Dispersion  Curves 


The  first  attempt  to  understand  the  short-period  wave 
trains  seen  in  Seneca  Lake  used  a  linear  theory.  A  matrix 
method  similar  to  that  of  Haskell  (1951)  was  developed  and 
a  Fortran  program  was  written  to  obtain  dispersion  curves 
for  internal  waves  in  a  many-layered  fluid.  As  work  progres¬ 
sed  on  the  Seneca  Lake  data  it  became  increasingly  evident 
that  the  wave  trains  seen  could  not  be  treated  as  linear 
waves,  so  this  approach  was  abandoned.  It  is  felt,  however, 
that  the  techniques  that  were  developed,  especially  the  For¬ 
tran  program,  can  be  of  use  in  other  problems  and  are  pre¬ 
sented  below. 

II. a.  To  derive  a  differential  equation  for  internal  waves 
The  derivation  of  the  differential  equation  used  to  get 
a  period  equation  can  be  found  in  a  number  of  textbooks 
(e.g.,  Phillips,  1966).  Starting  with  a  fluid  in  an  unper¬ 
turbed  state,  the  hydrostatic  equation  is: 


9Po 
3  z 


+  gp«  =  o 


Po  =  Po (z) 


Po  =  Po  (z) 


(II. a. 1) 


The  linearized,  frictionless  equations  of  motion  are  (neglec¬ 
ting  second  order  terms) : 
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3u 

3 1 

fv  = 

1 

Po 

9p 

3x 

(II . a. 2) 

9  v  , 

fu  = 

1 

3p 

3 1  + 

Po 

3y 

(II . a. 3) 

9w 

9 1 

-p-  P 

P  0 

1 

Po 

3p 

3  z 

(II. a. 4) 

with  f  =  2  ft  sin  X  (  ft  the  angular  velocity  of  rotation  of 
the  earth  and  X  the  latitude) . 

The  fluid  is  assumed  to  be  incompressible: 


DP  _  n_  3p  ,  3po  TT  (II.  a.  5) 

-  —  U  —  TT—  +  — W 

Dt  9t  9z 

The  equation  of  continuity  is: 


3u  +  _3v  +  _3w  _ 
3x  3y  3z 


(II. a. 6) 


When  the  fluid  is  in  a  perturbed  state  the  density  and  pres¬ 
sure  will  be: 


P 

ii 

ts 

o 

+ 

ts 

p' 

=  p' 

( t ,  x  ,y  , 

-  z) 

P 

CL 

+ 

o 

Cl 

II 

p' 

=  p' 

(tfX.y, 

-  z) 

Equations  (II. a. 2)  through  (II. a. 5)  become: 


3u 

3 1 

-  fv  =  - 

P  0 

3p ' 

3x 

(II. a. 2a) 

3  v 

3 1 

+  fu  =  - 

P  0 

3p 1 

3y 

(II .  a .  3a) 

9w 

3 1 

+  g-£-  =  -  — 
P  o  P  0 

3p ' 

3  z 

(II. a. 4a) 

3  P  ' 

3 1 

+  ^p0  w  -  0 

3  z 

(II. a. 5a) 

Equations  (II. a. 4a) 

and  II. a. 5a) 

are  used  to  eliminate  : 

Po 

9 

32w  3pn  1 

It7--  3  z  °  w-  g 

3  2p 7 

3  z  3t 

(II .a. 7) 
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By  combining  equations  (II. a. 2a)  with  (II. a. 7)  and 
(II. a. 3a)  with  (II. a. 7)  two  equations  can  be  obtained: 


3  3w 

3 12  3x 

3  3  u 

3 12  3  z 

+  f  32v 
3t3z 

+ 

3x 

=  0 

(II .a. 8) 

3  3w 

3t 2  3y 

3  3  v 

3t2  3  z 

£  32u 
3t  3z 

+ 

3y 

=  0 

(II. a. 9) 

with 

N2  =  - 

g  9p0 

Po  9z 

N  is 

a  function  of 

z  only. 

has 

the 

dimensions 

of  seconds- 

and  is  called  the  Brunt-Vaisala  frequency. 

With  the  equation  of  continuity,  both  u  and  v  can  be 


eliminated 

from 

equations 

(II 

.  a .  8) 

and  (II. a. 9) 

to  yield: 

3  f  32w 

32w 

3  2w  1 

2  ' 

r  32w 

3  2w  ~f 

3t 2  L  8x2  ' 

3Y  2 

+  3z2J  + 

$  i 

L  3x2 

'  3Y  2  J 

(II. a. 10) 

It  is  assumed  that  a  wave  solution  for  the  vertical  velocity 
w  in  the  horizontal  plane  exists: 


w(x,y,z,t)  =  e^  kxx  kyy^(z) 

with  k„  and  k  the  wave  numbers  in  the  x  and  y  directions, 

x  y 

Equation  (II. a. 10)  then  reduces  to: 


(m2-N2)  (k  2+k  2 )  cf)  +  (f2  -03  2  )il£  =  0 
x  Y  3  z 2 


(II .a. 11) 


Rearranging  terms  and  calling  the  wave  number  in  the  direc¬ 
tion  of  travel,  k,  the  differential  equation  is  obtained: 


Hi  ,  ,.2  u2-N2 


3  z' 


+  k 


f  2  — GO 


(p  =  0 


(II. a. 12) 


with  k  =kx  +ky 


Il.b.  To  obtain  a  period  equation  for  the  many-layered  model 
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using  a  matrix  method 

The  derivation  of  a  period  equation  for  internal  waves 
in  a  fluid  of  many  layers  using  a  matrix  method  is  similar 
to  a  method  developed  by  Haskell  (1951)  for  surface  seismic 
waves.  We  start  with  a  fluid  of  finite  depth  with  a  rigid 
bottom  and  take  the  z  axis  pointing  positively  down.  In  gen¬ 
eral,  fj  ,  the  Brunt- Vaisala  frequency  will  be  a  continuous 
function  of  the  depth,  z.  This  will  be  approximated  by  assu¬ 
ming  the  fluid  to  contain  n  layers,  each  with  a  constant  Vai¬ 
sala  frequency  and  a  discontinuity  of  the  Vaisala  frequency 
across  layers.  We  wish  to  solve  the  differential  equation: 


(Il.b.l) 


with  [J  constant. 


In  layer  p  we  have: 


0 


(II. b. 2) 


whose  solution  is: 


(II. b. 3) 


with 


x 


(II. b. 4) 


For  03>  f  ,  ir\p  will  be  real  for  Np  and  imaginary  for  wtfjp  . 


Taking  the  derivative  of  4>  with  respect  to  z,  which  will  be 


denoted  (j>,  we  obtain: 


(II. b. 5) 
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Layer  p  with  thickness  h  is  bounded  by  the  (p-1)  inter- 

P 

face  on  top  and  the  p  interface  on  the  bottom.  Coordinates 
can  be  translated  so  that  the  origin  of  the  z  axis  is  on  the 
(p-1)  interface  (this  only  changes  the  constants  Ap  and 
which  will  drop  out  of  the  results) .  At  the  (p-1)  interface 
we  have : 


^p-1  “  ^  +  Bp 


(II. b. 6) 


Vi  =  WVp 


(II. b. 7) 


This  can  be  written  in  matrix  form  as: 


( cb 

VHP"1, 


Vi 


>  = 


W 


V 


(II .b. 8) 


with  E  the  matrix: 

P 

- 

1  1 

E  = 

[  "V  _mp. 

Written  out,  equation  (II. b. 8)  is: 


7 

i 

r  n  1 

r  n 

*P-l 

1  1 

ap 

• 

^p-1 

m  -nr 

p  pJ 

bp. 

• 

At  the  p  interface,  z-h^  anc^: 

*p =  vmphp  +BPe_mphp 

»P 


or  in  matrix  form: 

( /  4*p)  —  Dp  ^ ^p '  ^p^ 


(II. b. 9) 


(II .b . 8a) 


(II. b. 10) 
(II .b. 11) 


(II. b. 12) 
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with  D  the  matrix: 
P 


dp  = 


eVp 

“HI  ll 

e  P  P 

m  h 

-m  h 

V  P  P 

"V  p 

( II  .b  .  1 3) 


Combining  equations  (II. b. 8)  and  (II.b.12)  to  eliminate 


A  and  B  results  in: 

ir  Jr 


_1  *  r  .  ( II . b . 14 ) 


(Ap,  Bp)  =  Ep  (  Vl'  Vl)=  DP  (V  V 


or 


^p'  ^p}  “  DpEp  ((})p-l'  ^p-l} 


(II. b. 15) 


with : 


-1 


1  1 
2  2 


mp 


1 

2 


2  m 


P 


We  define  a  new  quantity,  a^DpEp^: 


i  i 


2  1 


1  2 


2  2 


( II . b . 16 ) 


( II .b . 17) 


The  terms  of  the  matrix  ap  will  be  hyperbolic  sines  and 
cosines  when  mp  is  real  and  trigonometric  sines  and  cosines 
when  mp  is  imaginary.  As  an  illustration,  consider  the  matrix 


element  ( ap)  : 


/  \  1  ,  ^  t")  I  m-p*  h  -pv  \ 

(a  ) =  i(e  P  P  +  e  P  P) 


P'  i  i  2 
For  rn  real: 


(II. b. 18) 


(V.. 


cosh (mphp) 


( II .b . 19 ) 
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For  nv  imaginary: 


We  therefore  have  the  equation: 

(V  V  =  VVi'Vt 


(II  .b. 20) 


( II . b  .  2 1) 


(II.b.22) 


In  a  similar  manner  we  can  obtain  for  the  (p-1)  layer,  of 
thickness  hp_^  and  Vaisala  frequency  ^  and  bounded  by  the 

(p-2)  and  (p-1)  interfaces,  the  equation: 


(<J> 


P 


-1' 


i 


(II .b. 23) 


In  equation  (II.b.22),  cf)  ^  an<3  ^D-q  refer  to  the  values 
of  <*)  and  <j>  of  layer  p  at  interface  (p-1)  ,  while  in  equation 
(II.b.23)  they  refer  to  the  values  of  (p  and  (f>  of  layer  (p-1) 
at  interface  (p-1).  Because  we  want  (as  boundary  conditions) 
both  <J)  and  }  to  t>e  continuous  across  interfaces  between  layers, 
the  values  of  (p  and  <J> -i  are  the  same  for  equations  (II. 

b.22)  and  (II.b.23)  and  we  do  not  have  to  denote  which  layer 
they  refer  to,  as  long  as  we  use  their  values  at  the  interfaces 
only.  Equations  (II.b.22)  and  (II.b.23)  can  be  combined  to 
obtain : 


(4> 


p' 


d>  )  = 


apap-l 


((J)p-2  ' 


^p-2 


) 


(II. b. 24) 
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with : 


a  a  n 
p  p-i 


(a  )  (a  ) 

P  1  1  P  1  2 


(aD)  U  ) 

P  2  1  P  2  2 


(a  i)  ( a  ) 

P  1  1  1  P  1  1  2 


aD_])  ( ap-l) 

P  ±  2 1  r  1  22 


=b 


(II. b  .25) 


If  this  process  is  continued  upward  to  the  first  layer, 
bounded  by  the  0  and  1  interfaces,  and  downward  to  the  last 
(n)  layer,  bounded  by  the  (n-1)  and  n  interfaces,  the  resul¬ 
ting  equation  is: 

^  ^n  '  V  =  anan-l  *  *  *  ap+lapap-l  ‘  '  a2ai  ^  (II.b.26) 


Call  the  matrix  product  of  a^ ,  F,  i.e.: 


F  =  a  a  n  ...  a  a  = 
n  n-1  2  i 


i  i 


2  1 


1  2 


2  2 


(II.b.27) 


Using  the  matrix  F,  we  have: 


=  *o> 


(II.b.28) 


4>  and  <f>  are  the  values  of  4>  at  the  bottom  and  sur- 
n  o 

face,  respectively.  c|>  is  set  equal  to  zero  because  of  the 
rigid  bottom.  cj)  0  is  assumed  to  be  zero  because  we  are 
dealing  with  internal  waves,  whose  surface  motions  are  very 
small . 


=  0 


(0,  (j>n)  =  F(O,4>0) 


(II.b.29) 


Equation  (II.b.29)  can  be  written  out 


0  =  0  F„  +  F12  <» 


(II. b. 30) 
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3>n  =  0  F  2  i  +  F  2  2  io  ( II .  b  .  31 ) 

Either  of  the  two  equations  can  be  used  as  a  period 
equation.  Equation  (II.b.30)  is  chosen  because  it  yields  the 
simpler  one: 

F 

1  2 


0 


(II. b. 32) 
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II. c.  Fortran  program 


IMPLICIT  KEAL*8  ( A-h,0-Z ) , REAL*4  ( W - Z ) 

REAL*4  SOI  (  5,200,4)  , ACNG ( 20Q » 4 )  ,C,PFR,  PERM 

01  Mb  NS  I  ON  DENS ( 51 ) , ZEN ( 50  )  , THKFT ( 50 ) , ENPRMN ( 50 ) , DATE ( 3 ) 

COMMON  EN( 50)  *  THKNS ( 50 ) , ML  AY 
COMMON/ F I XUP/FCUR, RTS( 50) , ABRTt  50 ) 

DATA  PI ,GRAV/3. 1415926536  0  0,0.950  D  3/ 

C  READ  INPUT  DATA 

200  R  E AD ( 5  ,  1  )  I  END , NPT S , NMCDE , NL A Y , I f T 
C  I  END  =  NO'.  TIMtS  THRU  LOOP  IN  SLVPf  NPTS=NO.  OF  PTS.  TO  CALCULATE 

C  66TWEI N  GMEGA=FCOR  AND  ENMAX  NMOOE=NO.  OF  MODES  TO  BF  CALCULATED 

C  NLAY  =  NO.  UF  LAYcRS  I F  T  =  1  FOR  THKNS  IN  FT.,  0  IF  IN  C". 

1  FORMAT  (5I=>) 

READ (5, 2)  T I  , DKK , OMGST , AK  \ST , AKFSH, AKPMN, DENSE 
C  TF  =  TOLFR  ENCE  OF  DEL  IN  SLVPE  DKK  =  DEL  T  A  KAPPA  USED  IN  DEl’E 

C  OMGZ  =  LOWEST  OMEGA  (PROG  WILL  CALCULATE  A  /APPA  FOR  OMEGA  BETWEEN  DM ..  Z  a\'D 

C  MAX  VAISALA  FRFO  (FNMAX).'  IF  OMGZ  .LT.  FCOR  WlLi.  SET  OMGZ=FCGR). 

C  AKAST  =  ST  ART  ING  (LOWEST)  i  -  A  P  P  A  AK  F  SH  =  H  I  GHi,  S  T  KAPPA 

C  TO  TEST  WITH  OMGST  AK PMN= M I N I  MUM  CHANGE  IN  KAPPA  FOR  NEW  MODE 
C  DENSB=DEN$.  A1  BOTTOM  OF  LA  T  (NLAY)  LAYER 

2  FORMAT ( 7F10.5) 

R  F  AD ( 5 , 3 )  (DENS(I) , THKNS (  I  )  ,  1  =  1, NLAY) 

C  DENS ( I ) =1  ENS.  AT  TOP  OF  LAYER  THKNS ( I ) =TH ICKNESS  OF  LAYER,  FI  0  C' 

3  FORMAT ( oFlO.5) 

RE  AD ( 5 , 99 )  DATE 

C  DATE  IS  IN  FpRM  (  I ) = Y  R  (2)=M0  (3)=uAY  (4)=HR  *(5)=M1N 

99  FORMAT (5FIQ.3) 

C  LADITUDE  AND  CORIOLIS  PARAM. 

RE Au ( E , 3 )  ALTDG, ALTMN, ALNDG, ALNMN 
ALATDG=aLT0G+ALTMN/60.0 

FCOR= ( PI / ( 12.0  D  0*3.6  D  3  )  )  *CS I N ( AL ATDG*P I / l BO . 0  D  0)*2.0  D  0 

IF  (NLAY)  1000,1000,20 

20  NLAY1=NLAY+1 
DENS (NLAY1 ) =  D  E  N  S  3 
IF  (IFT)  24,24,22 

22  DO  23  1=1, NLAY 

23  THKNS  (  I  )  =  1?.0  D  0*2.54000^  D  (j*THKNS(I) 

C  CALCULATE  VAISALA  FREQS,  FIND  HIGHEST,  SET  SOLUTION  ARPAY  TO  ZbRO 

24  DO  124  1=1, NLAY 

124  L  N (  I  )=DSORT (GRAV*DLOG( DENS ( 1  +  1  )/DENS (I)  ) /THKNS ( I  )  ) 

CALL  MAX ( - M  M  4  X  » E  N , N  L  A  Y ) 

JN=NPTS+ 1 
DO  21  1=1,5 
00  21  J  =  l  ,  J N 
DO  2  1  K  = l , 4 

21  SUL ( I , J , K ) =0 . 0  D  0 

C  FIND  OMGST  (OMEGA  HALF  WAY  BETWEEN  OMGZ  AND  ENMmX  ON  LOG  SCAL.).  U3:S  THIS 
C  TO  FIND  THE  FIRST  KAPPA 
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IF  (  FT.  OR  .GT.  OMCZ)  OMGZ-FCOR 
DL,OMG=DABS  (OLOG(  ENMAX/OMGZ  )  ) 

DOMG=DLOMG/NPTS 

OMGL=DLOG( OMGZ ) +NPTS*D0MG/2.0  0  0 
OMGS  T=DEXP ( OMSL ) 

C  PRINT  AND  PUNCH  INPUT  DATA 
WHITE (6,4) 

4  FORMAT (  •  1 •  ,  1  5  X  , '**  INTERNAL  WAVE  DISPERSION  IN  FLUID  OF  MANY  LAYERS 
1*  EACH  WITH  CONSTANT  VMSALA  FREQUENCY**'  ,//,2f)X,  MJSFS  MATRIX  M  TH 
20D  SIMILAR  TO  THAT  OF  HASKFLL  TO  SET  UP  PERIOD  F QUA  I  I  ON • , / , 40 X , 

3  'SOLVES  BY  ITERATION  OF  KAPPA  FOR  A  GIVEN  OMEGA*,//) 

WRITE  (6,9.8  )  DATE,  ALT  DG,  AL  T MN , ALNOG , A LN"N 
98  FORMAT ( T45, * *DATE  AND  POSITION  OF  MOOFL* ' , / , T 10, ' YE AP ' , F6. C , 3X , 

I  'MONTH' ,F4.0,3X, 'DAY' ,F4.0,1UX, 'LATITUDE* ,F7.0, '  UEG',F4.0,'  MIN' 
6  , 10X, • LONGITUDE' , F7.0, •  DEG',F4.0f'  MIN',//) 

WRITE (6,5)  FCOR, OMGST,AK AST, N LAY, N MODE 

5  FORM AT ( / , 1  OX , 'CORIOLIS  PARAMETER  = ' , D 1 ? . 5 , 5X , ' STAR T I  MG  OMEGA  =', 

1  D12.5, 5X, • START ING  KAPPA  -• , 012*5, /,20X, 'NUMBER  OF  LAYERS  USED  IN 

2  MODEL  =  •',  13, 5X, 'NUMBER  OF  MODES  TO  BE  CALCULATED  =',I2,/) 

ZFCOR=  FCOR 

ZDNSB=0CNS8 

ZLAY=NLAY 

WR I TE ( 7,9  7  )  DATE , ALTDG, ALNDG, ZFCOR , ZDNSR, ZLAY 
97  FORMAT ( 5FS. 1, 3EI3.6) 

WRITF(6,6) 

6  FCRMAT(/,T2U, 'THICKNESS  OF  L A YFR • , T69 , ’ TOP  DENS.  OF  LAYER', T100, 

1  'LAYER  VAISALA' ,/,T2,'LAYEis  Nfl.',Tl9,'IN  CM.',T37,'IN  FT.',T73, 

2.  'IN  GM./c.C. ' ,T9I, ' FREQ.  IN  I/SEC .', T 108 ,' PER.  IN'  MIN.') 

DO  223  1  =  1*,  NL  A Y 

THKFT ( I )=1 HKNS( I )/( 12.0  0  0*2.540004  D  0) 

Z EN ( I )  =  E N (  I  ) 

223  ENPRMN(  I  )  =2.0  D  0*  P  I  /  (  EN  (  I  )  *60 . 0  C.  0) 

W  K I  I F ( 6 , 7  )  ( I ,THKNS(  I  ) ,THKF1 ( I ) ,DFNS( I ) ,EN( I ) , ENPrMN ( I ) , 1  =  1 ,NLAY) 

7  F  0  RMAT(T4, I3»T14,  Tj  13. 6, T32,L13.6,T  70,  17.10, T  9 1,D  13.6, T10B,  013.6) 

WRITE (6, 107)  DENS ( NL  AY  1 ) 

107  FORMAT ( T70,H17. 10) 

WR I TE (7,94  )  ( THKNS (  I  ) , DENS ( I ) , l EN ( I ) , I  =  1 ,NLA Y ) 

94  FORMAT (2(F12.2,F12.7,Ll6.tt) ) 

MODE= 1 
NPTM=0 

C  FIND  KAPPA  FOR  OMGST  MP=1  FOR  SOLUTIONS  .GT.  OMGST 

25  OMG=OMGST 
N  P  T  N  =  C 

MP  =  1 

CALL  ROOTS (NLA Y,OMG, FCOR, RTS, A BRTiEN) 

CALL  FRSTK(AKAST,AKPMN» AKAP,OMG,AkFSH,TF*DKK, I ER , I  END, F ) 

IF  (IER)  154,26,154 

154  WRITF(6,155)  IER 

155  FORMAT ( / , T  5  0 , '**EHROR  RETURN** • , / , T55 , • I EP  =  • , I  2 , /) 

GO  TO  54 

C  GET  SOLUTIONS  ALONG  DISPERSION  CURVE 

26  0MG1=0MG 
AKAST=AKAP 

27  AKAPl=AKAP 
F  1  =  F 

0MGL2  =  DL0G(  OMiGl  ) +D0MG*0.C5  D  0 
0MG2=DEXP ( 0KGL2 ) 
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IF  (0MG2-ENM AX  )  127,30,30 
127  CALL  ROOTS  INLAY, 0MG2,FC0R, RTS, ABRT,EN) 

C-ALL  SLVPF  (0MG2,AKAP,TF,DKK,  ILR,  IEND,F  ) 

IF  (IER)  37,2°, 37 
28  AKAP2=AKAP 
F2  =  F 

0MG=(OMGH-OVG2)/2.O  0  0 
AKAP=(  AKA=>l  +  AKAP2  )/2.0  I)  0 
U= (0VG1-0MG2 )/<AKAPl-AKAP2) 

NPTN=NPTN+ 1 

IF  ( NPTM  .LT.  NPTM )  NPTM=NPTN 

SOL  (  MODE, iNIPTN,  1  )=0MG 

SOL (M0DE,NPTN,2 )=AKAP 

SOLI  MOD Ft  NPTN,3)=U 

SOL  I  MODE  » NPTN » 4 )  =  ( F 1  +  F2 ) / 2 .0  0  0 

OMGNL=DLOG(OMG)+MP*DONG 

0MG1=DEXP(0MGNL ) 

C  IF  OMEG A=  LNM AX ,  FINO  SOLUTIONS  .LT.  OMGST 
IF  I OMG l- P  NO AX )  32,30,30 

C  MP=- 1  FOR  SOLUTIONS  .LT'.  OMGST 
30  M P  =  -  1 

NPT  A  =  NPTN 
OMG=ONGST 

OMGNL=DLOG(OMG)+Mp*DOMG 
OMGI=DEXP (OMGNL ) 

U=  SOL  I  MODE ,1,3) 

AKAP=AKAST 

C  U-SE  SLOPE  OF  DISPERSION  CURVE  TO  GUESS  KAPPA  OF  *  NE  X  f  OMEGA 
32  IF  (0MG1-FC0R)  38,38,34 
34  AKAP=AKAP+(0MG1-0MG)/U 

CALL  ROOTS  INLAY, OVG l, F COR, RTS, APR T,EN) 

CALL  SLVPF (0MG1 , AKAP,TF,DKK, I ER , I  END , F  ) 

IF  (IER)  37,27,37 

37  IF  (MP)  38,84,30 

C  REARRANGE  SOLUTION  IN  ORDER  OF  INCREASING  OMEGA 

38  NPT6=NPTN-NP r A 
DO  4  2  1  =  1,  NPT *3 
J=NPTN- I + 1 

DO  42  K=l,4 

42  AKNGI I ,K)  =  SCL(MODE,J,K  ) 

DO  44  I  =  1 , NP  T  A 
J=NPT8+I 
DO  44  K=  1,4 

44  AKNGI J,K)=SOL(MODE,I,K) 

DO  46  1=1, NPTM 
DO  46  J=  1 ,4 

46  SOL  I  MODE  ,  1  ,  J  )  =  ARNE- 1  I  ,  J  ) 

C  PRINT  SOLUTION 

WRITE (6,8)  MODE 

8  FORMAT!/////, SOX, •*  MODE  NUMB t R ',13,'  *»,//, T6 ,' OMEGA ', T 1 V ,  • KAPP A ' 

1  ,T29,  ‘GROUP  VhL  U*  ,  T42,  •  PHASE  VcL  C  '  ,  T  5  7  ,  •  PE  R  1  OD  •  ,  T  70  ,  '  P  F  R  I  OP  •  ,  T  J l 

2  ,  ' WAVE  LENGTH* , T94, • WAVE  LENGTH* , Till, *F12,,/,T3,  • ( 1/SFU)  •  ,Tlv, 

3  *  I  1/CM  )  •  ,  T3U,,*  I  CM/ SEC  )  *  ,  143,  •  (CM/SEC-1  *  ,T53,  »  I  SEC  )  •  ,  T71  ,  *  I  MI  U)  •  , 

4  T  8  '  (CM)  •  ,  T98  ,  ’  (  FT  )  •  ,/  ) 

DU  80  I  =  1  , NpTN 

C=  SOL  I  MODE,  I  ,  D/SOLIMHCE,  1,2)  * 

PLR  =  2.0  l'  0*P I /SOL  I  MODE,  1 ,  1 ) 
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PFP' =PER/bO.O  D  0 
W  L  =  2  .  C  0  0*P I / SOL ( MODE  » 1,2) 

WLFT=WL/( 2. 540005  D  0*12.0  l  u) 

50  WR  I  TP (6, 10)  ( SOL ( I  ODE, I , J ) , J  =  1 , 3 ) ,C, PFR, °FRN  ,V.L,WLFT,SOL(  VObE,  1,4) 

10  FDR  NAT  (vJE13.r>) 

C  TEST  IF  MORE  MQuES  WANTED 
IF  ( *ODE-,NMOOE  )  52,54,54 

52  MODt=POOF+l 

AKPMN=AKPMN+  AK  AST 
GO  TO  25 

C  PUNCH  SOLUTION  ON  CARDS 
54  N P  T  S  S  =  N  P T  ’■  +  1 

Wr'  I  1  E  (7,93)  (  (  SCL  (  I  ,  J  ,K  )  ,K=1  ,  3  )  ,  I  =  1  ,  NMOD  F  )  ,  J  =  1  ,  N'P  T  S  S  ) 

93  FORfAT (3E13.o,2X»3E13.6) 

GO  TO  200 
1000  CALL  EXIT 
FND 


SUBROUTINE  FH.  STK ( AKAST ,AKPMN,AKAP,ONG,AKFSH,TF,LKK , IcR, IEND,F) 

C  FINDS  SMALLEST  KAPPA  G  - .  E  A  T  cR  TFiAM  THAT  SPECIFI  D  THAT  SOLVES  PERIOD  rQ. 

IMPLICIT  E  A  L  *  3  (A-H.P-7) 

1  AKAPs AKAST 

CALL  SLVPE(0MG,  AKAP  ,  TF  ,  DKK  ,  I  E  V  ,  I  f.ND,  F  ) 

IF  (  ICR  )  4,2,4 

2  IF  (AKAP-AKPKN)  a, 4, 8 

4  A*  A  S  T=  5 • 0  D  0  *  A K  A  S  T 

IF  (  AKAi>T-AKFSH)  1,1,17 

5  AKAST=AKPMN+0.2  D  0 * ( A K A S T -AK PMN ) 

DU  15  1=1,5 

A  K  A  P  = I *AKAST 

CALL  SLVPt  (CJMG,  AKAP,TF,DKK,  IER,  IEr'D,F) 

IF  (IER)  15,10,15 
10  IF  ( AK AP-AKPMN)  13,15,13 
15  CONTINUE 

17  I  F  R  =  2 

18  RETURN 
FND 


SUP.  TOUT  I  N r  v  AX  (  A «  A X  ,  A  ,  N  ) 
REAL*0  A  M  A  X , p , A (  1  ) 

A  iV  A  X  =  A  (  1  ) 

DO  10  1=2,3 
P  =  A  (  I  ) 

IF  (  A  k  A  X  .LI.  3)  A  Ni  A  X  =  H 
10  CONTINUE 
P  F  T  U  -  i 
END 
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SUBROUTINE  SIVPfc (  CMG »  AKAP, TF , DKK  , I  ER  ,  I  END , OE L ) 

C  FINDS  A  KAPPA  FOR  GIVEN  OMEGA  THAT  SOLVES  PERIOD  EQ.  USES  OFLPE  TO  TEST 
C  KAPPA 

IMPLICIT  REAL *8  (A-H.O-Z) 

I  ER  =  0 

HDKK=DKK 

TOL  =  TF  / 1.0  Lj  3 

CALL  DELPE(OMG*AKAP,DC-L,DF,DKKt  IFK) 

IF  (  IER  .GT.  OF  GO  TO  20 
DO  10  I  =  1 » I  END 
IF  <DF)  St 2,8 
2  D K K  =  2 • 0  D  0  *UKK 
GO  TO  35. 

5  DK  =  2 . 0  D  0  *DKK*AKAP*DEL/DF 

IF  (  DABS(  DEL  )  -DAPS  (  OF  /  2  •  0  l*  0))  12,15,15 
12  DKK=DKK/5.0  D  0 

IF  (DKK-0.5  D-15)  40,32,32 
15  DK  =  0 . 8  D  0*DK 

DKA  =  DA BS  (  0 •  5  D  0*>>K/AKAP) 

IF  (OKA  .LT.  DKK)  PKK=DKA 
IF  (DKK-0.15  D-14)  40,22,22 
22  IF  (DABSHDK/AKAPJ-l.O  D  0  )  32,25,25 
25  IF  (UK)  30,30,27 
27  AKAP  =  AKAP/3.0  D  0 
GO  TO  35 

30  AKAP=2.0  0  0*AKAP 
GO  TO  35 
32  AKAP= AKAP-DK 

35  CALL  DELPE (OMG, AKAP, DEL, DF,UKK, IER  ) 

IF  (IER  .GT.  OE  GO  TO  20 
TOLF=TOL*AKAP 

IF  (DKK-1.0  D-10 )  38,38,10 
38  IF  (DABS(DFL)-TOLF )  20,20,10 
10  CONTINUE 

C  ERROR  RETURN  NO  SOLUTION  WITHIN  TOLERENCES  AFTER  I  END  TIMES  THRU  LOOP 
IEP=1 
GO  TO  20 
40  TOLR=TF*AKAP 

IF  ( DABS ( DEL )-TOLR)  20,20,48 
C  ERROR  RETURN  SOLUTION  NOT  WITHIN  TOLERANCES 
45  I E  R=  3 
20  DKK=HDKK 
RETURN 
END 
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SUBROUTINE  DELPE ( OMG,  AK  AP , DEL  »  DF  »  nKK  »  IER  ) 

FINDS  DELTA  FOR  PERIOD  EO.  F( 1 , 2  )  =DEL  T  A  *  WHICH  WE  WANT  =  0 
F=A (N) AIN-l  )  . . .A(  I ) ,  A=  2  BY  2  MATRIX 
IMPLICIT  PEAL*8  IA-H,0-Z) 

DIMENSION  A(2,2),6(2,2),F(2,2),F1?(2> 

COMMON  EN(50) ,D(50),NLAY 

COMMON/  FIXUP/FC0R»RTS(50)  ,  AB.RTI50) 

ROOT  =  R  T  S ( NL  AY ) 

S<?0UT=ABRT(NLAY> 
n J=D ( NL AY ) 

AKAP=AKAP*( 1,0  D  O-DKK ) 

DO  100  11=1,2 

CALL  AMTRX (ROOT,AKAP,DJ,F,SROOT, IER) 

N=NL AY-1 
DO  15  1  =  1, N 
J=NLAY-I 
0 1  =D ( J ) 

RJOT=RTS ( J ) 

SRJ0T=A3RT ( J ) 

CALL  AMTRX(RJOT,AKAP,UI,A,SkJOT, IER) 

DO  10  K= 1 , 2 
DO  10  L=  1 , 2 
10  P(K,L)=F(K*L) 

DO  15  K= 1 , 2 
DO  15  U=l,2 

15  F(K,L)=B(K,1)*A(1,L)+.(K,2)*A(2,L) 

F121II  )  =  F(  1,2) 

100  AKAP=AKAP* ( 1.0  D  O+DKK ) / < 1 . 0  D  O-DKK) 

DEL=(E12(l)+F12(2))/2.0  0  0 

AKAP  =  AKAP*(  1.0  D  O-DKK  )  /  (  (  1  •  0  D  O+OKK  )  *  (  1  •  0  D  O+OKK )  ) 

DF=< F12(2)-F12( 1  )  ) 

RETURN 

END 


SUBR  OUT I NH  ROOTS ( NLAY, OMG , FCOK,RTS,ABRT , EN ) 

C  CALCULATES  ROOTS  IN  EACH  LAYER 

REAL*8  OMG,FCOR,RTS(  1  )  ,  A B k T  (  1  )  -. E N  (  1  ) 

DO  10  1=1, NLAY 

R  T  S ( I )= (OMG*OMG-FN( I  )*EN(  I  ) ) / (0MG*0*G-FC0K*FC0-  ) 
10  ARRTU  )  =  DSQ'RT  ( DABS ( RTS (  I  )  )  ) 

return 

END 


140 


SUBR0UT  I  NE  AMTRX  (  ROUT  ,  AKAR  ,  D,  A,  SHOOT,  I  EH  ) 

C  FINDS  MATRIX  ELEMENTS  FOR  INTERNAL  WAVES  IN  LAYERED  FLUID  USING  HA,SKcl 

C  METHOD 

REAL*'!  ROGT,AKAP,Uf  S  ROOT  ,  EM  ,  EM  D,  ET  ,  A  (  2,2) 

EM= AK \P*S°OOT 
EMD=EM*D 

IF  (ROOD  10,10,15 
10  A( 1,1 )=DCOS(EMD) 

A (  1, 2)=DSIN(E"D)/FM 
A(2, 1  )=-EM*EM*A( 1,2) 

GO  TO  20 

15  IF  (EMD  .GT.  1.74  D  2)  GO  TO  25 
ET  =  OEXP(Ei-'D) 

A(1,1)=(ET+1.0  D  0/ET)/2.0  D  0 
A-(l,2)  =  (FT,r.1.0  D  0/ E  T ) / ( 2  •  0  D  0*EM) 

A ( 2 , 1 )=EM*EM*A( 1,2) 

20  A ( 2 , 2 )  =  A ( 1,1) 

RE  TURN 
25  I  E R  =  5 
RETURN 
END 
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